Keywords: Antarctic fur seal, Arctocephalus gazella, abundance, assessment, decline, reduction, IUCN, Red List.


1 Introduction

The Antarctic fur seal (Arctocephalus gazella) is one of the smallest otariids (eared seals), and has a circumpolar distribution around Subantarctic and Antarctic islands of the Southern Ocean where it aggregates to breed. It has a four-month lactation, which is one of the shortest in its family, probably as an adaptation to the extreme seasonality of its nursing environment. Contrary to other Antarctic pinnipeds, which primarily rely on stored energy reserves during lactation, it requires constant foraging to rear the offspring, and heavily relies on a nearby productive and predictable foraging environment. Additionally, it is highly sexually dimorphic, with a polygynous breeding system that significantly determines longevity and lifetime breeding success. Despite being physiologically reproductive at ages three to five, mature males mostly hold breeding territories between ages nine and ten and have brief reproductive lives subsequently. Starting at age three and above, females normally produce a single pup nearly every year and can live up to 20 years or more. Due to its high fecundity and distribution in highly productive waters, the Antarctic fur seal is the second most abundant pinniped in the Southern Hemisphere, and the most abundant otariid globally (Wickens and York 1997; Costa, Weise, and Arnould 2007; Forcada et al. 2023b). As a numerically important predator of Antarctic krill (Euphausia superba), it is a main component of Southern Ocean ecosystems, particularly of the South Atlantic region (Croxall and Prince 1979; Croxall et al. 1988; Agnew 1997), being a key indicator species of their state, function, and services.

According to Humble et al. (2018) and Paijmans et al. (2020), there are at least four distinct genetic groups or subpopulations, with the largest found in the South Atlantic sector of the Southern Ocean, including South Georgia, Bouvetøya, and the South Shetland Islands (Forcada et al. 2023b; Forcada and Staniland 2018) (Fig1). Together, these three subpopulations comprise approximately 98% of the global abundance, but all have experienced rapid declines due to increasing external pressures (Forcada 2021; Forcada et al. 2023b; Forcada and Hoffman 2014; Krause et al. 2022, 2024; Hofmeyr 2016).

At Bird Island (NW South Georgia; Fig.2b), fur seal pup production declined from 60,500 to 27,900 between 2009 to 2017, or 9.2% per year, and the mature female population declined annually by 7.2% (Forcada et al. 2023b). At Bouvetøya (Fig1; 5), the second largest subpopulation with approximately 66,000 seals by 2002 (Hofmeyr et al. 2005), pup counts declined annually by 5.6%, between 2001 and 2007 (Hofmeyr 2016). And at the South Shetland Islands (Fig1, 4; Fig2a), with approximately 41,000 seals in 2002 (Krause et al. 2022), annual pup counts at study areas in the southwest declined by 90.9% between 2008 and 2023 (Krause et al. 2024), or a 5.6% per year. At the entire archipelago, counts declined from 10,057 pups in 2002 to 7,602 in 2008 (Krause et al. 2022), equivalent to an annual decline of 4.6%.

In the South Atlantic region, subpopulations are declining because of directional environmental effects, particularly at South Georgia (Forcada, Trathan, and Murphy 2008; Forcada and Hoffman 2014; Forcada et al. 2023b). However, there are known competing effects such as local top-down control (Schwarz et al. 2013; Krause et al. 2024), and other potential effects which are difficult to quantify (Forcada 2021; Forcada et al. 2023b; Krause et al. 2024). These include, for example, competitive exclusion from staple prey by increasing whale populations, e.g., Baines et al. (2022), Biuw et al. (2024), and Savoca et al. (2024), and by locally increasing commercial fisheries (Lowther et al. 2020; Watters, Hinke, and Reiss 2020); emerging pollutants (e.g., Aznar-Alemany et al. 2019), including microplastics (Garcia-Garin et al. 2020), with concerns over their persistence and toxicity (e.g., Brault et al. 2013); and outbreaks of new epizootic diseases such as avian influenza (H5N1) (Bennison et al. 2024; Brownell et al. 2024).

Other breeding locations, currently part of a single genetic group named here the ‘East Antarctica’ subpopulation, are distributed in subantarctic islands of the southern Indian and southwestern Pacific Oceans (Fig1; locations 6-10). Recently, pup counts at most locations, including Marion Island (Prince Edward Islands) (Wege et al. 2016), Possession Island (Crozet archipelago), Courbet Peninsula (Kerguelen archipelago) (Guinet et al. 2025), Heard Island (Heard Island and McDonald Islands) (Page et al. 2003), and Macquarie Island (Goldsworthy et al. 2009) (Fig. 1), have remained stable or slowly increasing annually, by approximately 1 to 4%.

The latest assessment of the species by the International Union for the Conservation of Nature (IUCN) (Hofmeyr 2016) recognises that the global population is in decline, which is mostly driven by the rapid declines observed in the subpopulations of the South Atlantic region. However, the Red List category assigned is Least Concern, which is unlikely to reflect the real magnitude of recent population reductions. Here, I reanalyse the change of each subpopulation using the latest available data to provide an up to date assessment for the species. I first reanalyse data from Bird Island, South Georgia, as main source of correction factors and adjustments across the species range, because comparable data from other regions are mostly lacking. The results include adjustments to scale up pup production counts to total mature population; and estimates of generation time (generation length, hereafter), required to evaluate population reduction. The analysis also evaluates and projects the overall population change at South Georgia with the most recently collected survey data at selected locations. I also analyse count data from Signy Island, representative of the seasonal male migration from South Georgia to the South Orkney Islands, to track and confirm by proxy nearly five decades of population change at South Georgia. Next, I use the most recent population data available from all the other species breeding locations to evaluate conservation status according to IUCN categories and criteria (IUCN 2024).

2 Materials and methods

2.1 Study populations, distribution and abundance

South Georgia, in the Southwest Atlantic region of the Southern Ocean (Fig1), is by far the largest Antarctic fur seal subpopulation with an estimated 97-98% of the global pup production and 3.5 [3.1, 3.9] million seals in 2009 (Forcada et al. 2023b). The subpopulation includes small numbers of breeders in the South Orkney and South Sandwich Islands archipelagos (Fig.1), also located in the Scotia Sea (Fig.2a), although the genetic identity of the South Sandwich seals has not been assessed. The South Orkney Island group (Fig.2a) has a population almost exclusively composed of non-resident juvenile and adult males in post-breeding temporary migration, largely from South Georgia (Boyd et al. 1998; Kightley and Caldwell 1982; K. A. Jones et al. 2021; Laws 1973), although there are small numbers of temporary migrants from the South Shetland Islands (Lowther et al. 2020; Ouled-Cheikh et al. 2024), given the relative proximity of both archipelagos. Males from South Georgia and South Shetland Islands typically migrate temporarily along the Weddell Sea and west Antarctic Peninsula, and as far South as the Bellingshausen Sea (K. A. Jones et al. 2021; March et al. 2021). There are also incidental records of pups in the Antarctic Peninsula, e.g., near Palmer Station (Larsen et al. 2025; Parmelee et al. 1977), but these are exceptional and extralimital within the species current breeding range.

The numbers of breeding females and pups at the South Orkney Islands are extremely low, and counted as less than 70 in the entire archipelago in the early 1970s. At Signy Island, which has been monitored annually since 1977 (Kightley and Caldwell 1982), despite several annual counts of more than 20,000 individuals between 1994 to 2000, only a maximum of five pups were recorded in any year of the same period, and pup numbers have declined to near zero since 2000 (Waluda, Gregory, and Dunn 2010). Similarly, at Laurie Island, only four pups were observed from 1996 to 2005, despite a maximum of 16,610 individuals counted in early 2005 (Carlini et al. 2006). At the South Sandwich Islands, breeders and pups have rarely been counted due to difficult access to locations, but since the early 1960s (Holdgate 1963) and after the regional recovery from sealing, pup numbers increased to at least 1,750 by 1997/98, of which 500 were attributed to Zavodovski Island (Convey, Morton, and Poncet 1999), which was counted again in 2023/24 (Table 1). Published records with more recent data for other islands, including Visokoi which had a count of 900 in 1997/98, are lacking.


**Figure 1**: Distribution of breeding populations of Antarctic fur seals with detail on genetic subpopulations. The table shows estimated pup production up to 2010 or before, for remote locations not recently surveyed. Colours corresponds to identified genetic subpopulations

Figure 1: Distribution of breeding populations of Antarctic fur seals with detail on genetic subpopulations. The table shows estimated pup production up to 2010 or before, for remote locations not recently surveyed. Colours corresponds to identified genetic subpopulations

**Figure 2**: **a** Scotia Sea, with breeding locations of the putative genetic subpopulation of South Georgia, including South Sandwich and South Orkney Islands, and the subpopulation of the South Shetland Islands. **b**, Locations of South Geogia where population surveys were obtained in 2021-22. In locations 2-5, seal counts were obtained with high resolution images from Remotely Piloted Aerial Systems (**RPAS**)

Figure 2: a Scotia Sea, with breeding locations of the putative genetic subpopulation of South Georgia, including South Sandwich and South Orkney Islands, and the subpopulation of the South Shetland Islands. b, Locations of South Geogia where population surveys were obtained in 2021-22. In locations 2-5, seal counts were obtained with high resolution images from Remotely Piloted Aerial Systems (RPAS)

2.2 Data

2.2.1 South Georgia

Antarctic fur seal demographic, genetic, and count data representative of South Georgia were collected at a long-term monitoring location of Bird Island named Special Study Beach (SSB) (54.0116\(^\circ\)S, 38.0507\(^\circ\)W; location 1, Fig.2b), from 1979 to 2025 (Croxall and Prince 1979; Doidge, Croxall, and Baker 1984; Duck 1990; Forcada, Trathan, and Murphy 2008; Forcada et al. 2023b). Additional data were obtained from Maiviken (54.2435\(^\circ\)S, 36.5033\(^\circ\)W; location 4, Fig. 2b), where a set of colonies named Burnet 1-3, Little, Poa 1-3 and Tortula are counted annually during the breeding season. Maiviken is situated on South Georgia main land approximately 105 km to the east of Bird Island.

2.2.1.1 Bird island capture-recapture and productivity data

Capture-mark-recapture of individuals, monitoring of pup production, and collection of individual count data at Bird Island used in this analysis started in 1979 (Doidge, Croxall, and Baker 1984) and have been consistently collected since 1984 onwards (Duck 1990; Forcada, Trathan, and Murphy 2008). Annually at SSB, from November 1st to mid-January, daily numbers of new and total females, newborn pups and territorial males are consistently recorded during twice-daily surveys. Within a season, female breeding status and performance is monitored through tags, and territory tenure by males is monitored through individual temporary paint marks. For between-year recaptures, a variable sample of females has been captured since 1984 every season, reliably aged from sectioned post-canine teeth, marked with plastic tags on the fore flippers, and subcutaneous electronic chips (Passive Integrated Transponder tags, PIT-tags) fitted since 1997. Flipper tags have been given to each cohort of newborn pups since 1984 and PIT-tags implanted since 2001; and genetic samples from pups and breeding females are collected upon capture. For each territorial male, genetic tagging is obtained each season with a small biopsy remotely acquired using a low-power crossbow bolt with a modified sampling syringe. Females and males are highly philopatric and return ashore to breed almost every year (Lunn and Boyd 1991; Hoffman and Forcada 2012; Hoffman, Trathan, and Amos 2006), once established as breeders and territorials respectively. Early pup mortality is estimated from the recovery of dead pups, prior to flipper-tagging approximately 1.5 months after birth.

2.2.1.2 Count data

Count data were collected as close as possible to the peak pupping day. The pupping season, from mid-November to mid-January, peaks with high spatial synchrony between the 6th to 10th of December, when numbers of newborn pups, breeding females, and territorial males at each location reach their maximum daily values, and show lower day-to-day variability than at later dates (Duck 1990; Forcada et al. 2005; Hofmeyr et al. 2007; Goebel et al. 2009).

Apart from the daily counts at SSB during the breeding season, complete counts of the entire Bird Island were obtained in nine seasons, including visual counts in December 1988, 1989, 1990, 1994, and 1997, and photographic counts in December of 2003, 2008, 2016, and 2023. (Hereafter, a breeding season is referred to by the second year of the split-year austral summer; e.g., 2023-24 is 2024.). Counts were conducted from elevated vantage points in all the breeding beaches by trained biologists.

Comprehensive aerial photographic data of South Georgia mainland were obtained during two helicopter surveys of all the fur seal breeding beaches conducted in December 2006 and 2008. For this analysis I use the count data available at the repository Forcada et al. (2023b), with emphasis on counts corresponding to locations 2 to 5 in Fig.2b to reassess the population at South Georgia. Surveys of these locations were repeated in December 2021 using RPAS (Remotely Piloted Aerial Systems). In these and other image-based surveys, two observers working independently using consistent training and methodology produced counts from digital images, classifying seals in four categories: pups, females, territorial males and non-territorial males.

From 2009 to 2025, colonies at Maiviken were visited every second day during the summer breeding season (November-February) to count numbers of fur seal males, females, pups and juveniles (Hollyman 2021). The data record is from 2009 to date (2025), although counts from 2009 were sparse and likely to be inconsistent in methodology compared with subsequent records. Therefore, for this analysis I only retain the data for 2010 to 2025.

All together, and with a maximum straight-line distance between any two counting or monitoring locations of 125km, the data obtained for 2021/22 provide a representative sample of geographic population variation over space and time at South Georgia.

2.2.2 South Orkney Islands

2.2.2.1 Count data

At Signy Island, South Orkney Islands (60.3385\(^\circ\)S, 45.0037\(^\circ\)W; Fig.2a), count data were collected around the entire island as part of a long-term monitoring programme which started in 1977 (Kightley and Caldwell 1982), and continued to date with only a few interruptions due to weather or logistic constraints (Waluda, Gregory, and Dunn 2010; Dunn et al. 2025). Counts by two trained observers using consistent methodology took place between mid-February and early March every year, when fur seals numbers peak.

2.2.3 Pup count data for other subpopulations

Available pup production and abundance data up to 2004 for all the species range are summarised in SCAR (2006); although some of the numbers reported are known to be unreliable because are only guesses, due to lack of proper counting efforts, and use biased adjustments, all of which needs revision and reassessment, as is the case of the subpopulation of South Georgia. For other regions, many reported records are total pup counts at breeding locations, from which extrapolations to total population size were also produced in some occasions. As discussed in the Population-reduction section below, earlier extrapolations are based on inaccurate scaling factors, and, also, do not separate mature population from pre-breeders, which is a requirement in Red List assessments. Here, I only provide total pup production counts or estimates, summarised in Table 1, which I used in subsequent calculations together with the new mature population estimates for South Georgia.

Breeding location season count method source Notes
South Sandwich Islands 1997/98 1,750 Visual Ref. 1 Late February; UC
Zavodovski Island, SSI 1997/98 500 Visual Ref. 1 Late February; UC
2023/24 1,265 RPAS Dickens (pc) Peak pupping
South Shetland Islands 1995/96 9,530 Visual Ref. 2 Late January; UC
2001/02 10,057 Visual Ref. 3 Late January; UC
2007/08 7,602 Visual Ref. 4 Late December; UC
2024/25 901 Visual Ref. 5 Late December; UC
Nyrøysa, Bouvetøya 1996/72 15,489 (754) Visual Ref. 6 Early January
1998/99 15,448 (370) Visual Ref. 6 Early January
2000/01 15,215 (353) Visual Ref. 6 Early January
2001/02 15,523 (480) Visual Ref. 6 Early January
2007/08 11,000 ? Ref. 7 Extrapolation?
2014/15 8,080 (310) Visual Lowther (pc) Early January
Marion Island 1999/00 343 Visual Ref. 8 February; UC
2000/01 464 Visual Ref. 8 February; UC
2003/04 744 (5) Visual Ref. 8 February; UC
2006/07 1,105 (9) Visual Ref. 9 January; UC
2009/10 1,379 (27) Visual Ref. 9 January; UC
2012/13 1,553 (64) Visual Ref. 9 January; UC
2015/16 2,052 (107) Visual de Bruyn (pc) January
2018/19 2,434 (169) Visual de Bruyn (pc) January
2021/22 2,528 (134) Visual Jordaan (pc) January
2024/25 3,172 (156) Visual Jordaan (pc) January
Prince Edward Island 2000/01 404 Visual Ref. 9 Late December; UC
Crozet, Possession Is 1999/00 234 Visual Ref. 10
2003/04 295 Visual Ref. 10
2020/21 504 Visual Ref. 11
2023/24 530 Visual Ref. 11
Kerguelen, Nuageuses Is. 2000/01 3,600 Visual Ref. 12 December; UC
Kerguelen, Courbet Peninsula 2000/01 1,600 Visual Ref. 12 December; UC
2009/10 4,287 Visual Ref. 11 December; UC
2023/24 4,793 Visual Ref. 11 December; UC
Heard Island 2000/01 1,012 Visual Ref. 13 November - March
2003/04 1,278 Visual Goldsworthy (pc)
McDonald Island 1979/80 100 Visual Ref. 10
Macquarie Island 1999/00 152 Visual Ref. 10 November - March
2003/04 165 Visual Ref. 10 November - March
2011/12 189 genetic Goldsworthy (pc) November - March

Table 1. Antarctic fur seal pup counts obtained at breeding locations other than South Georgia and South Orkney Islands. Parentheses in count are estimated standard errors, normally obtained from mark-recapture estimates of pup production. method is either visual, RPAS (Remotely Piloted Air System), or both. source is either literature cited, or personal communications (pc) from various contributors, including those listed in next section. References cited are as listed here; Ref. 1: Convey, Morton, and Poncet (1999); Ref. 2: Meyer, Walker, and Holt (1996); Ref. 3: Goebel et al. (2003); Ref. 4: Goebel et al. (2009); Ref. 5: Krause, Woodman, and Goebel (2025); Ref. 6: Hofmeyr et al. (2005); Ref. 7: Hofmeyr (2016); Ref. 8: Hofmeyr et al. (2006); Ref. 9: Bester, Ryan, and M. (2003); Ref. 10: SCAR (2006); Ref. 11: Guinet et al. (2025); Ref. 12: Lea et al. (2008); Ref. 13: Page et al. (2003). In Notes, UC stands for undercount.

2.3 Analysis

2.3.1 Population modelling

2.3.1.1 South Georgia

I obtained an up-to-date assessment for Bird Island and other locations at South Georgia with the latest available data and using population analysis methods of Forcada et al. (2023b). The analysis used integrated population models (IPMs) to evaluate population numbers, vital rates and other demographic parameters, and corrections for seal survey counts. It combined different data sources including individual recapture histories of fur seal pups and adults, pup production, and individual daily count data. To evaluate abundance change and trends at locations away from SSB, I assumed that the dynamics of their daily counts are highly correlated to those at SSB. This is accurate within Bird Island and with Maiviken (MVK), another study site located 105km East of Bird island as the crow flies (Forcada et al. 2023b), which is closer to the additional survey locations in 2022 (Fig.2b). However, I relaxed this assumption by also incorporating daily counts at Maiviken in the analysis. This accounted for differences in haul out patterns between this study location and SSB, and corrected counts at other locations based on their proximity to either study site.

2.3.1.1.1 Integrated population models (IPMs)
2.3.1.1.1.1 Females

I used the same models as described in Forcada et al. (2023b) with different formulations for females and males. For females, the expected number of individuals in year \(t\)+1 is described with the following age-stage-structured population projection model, which assumes a birth pulse and a pre-breeding census:

\(E([N_{t+1}|N_{t}])=E\left( \begin{bmatrix} N_{1,t+1}|N_{1,t}\\N_{2,t+1}|N_{2,t}\\N_{3,t+1}|N_{3,t}\\N_{4,t+1}|N_{4,t}\\ \vdots\\ N_{7,t+1}|N_{7,t}\\ N_{7+,t+1}|N_{7+,t}\\N_{N,t+1}|N_{N,t}\\N_{B,t+1}|N_{B,t}\end{bmatrix}\right)=\)


\(\begin{bmatrix} 0 & 0 &\varsigma_{F}\phi^0\phi^1\phi^3\alpha^3 &\dots &\varsigma_{F}\phi^0\phi^1\phi^7\alpha^7 &0 &\varsigma_{F}\phi^0\phi^1\phi^N\psi^{NB} &\varsigma_{F}\phi^0\phi^1\phi^B\psi^{BB}\\ \phi^1 & 0 & 0 &\dots &0 &0 &0 &0 \\ 0 & \phi^2 & 0 &\dots &0 &0 &0 &0 \\ 0 & 0 & \phi^3(1-\alpha^3) &\dots &0 &0 &0 &0 \\ 0 & 0 & 0 &\ddots &0 &0 &0 &0 \\ 0 & 0 & 0 &\dots &\phi^7(1-\alpha^7) &0 &0 &0 \\ 0 & 0 &0 &\dots &0 &\phi^{7+} &0 &0 \\0 & 0 &0 &\dots &0 &0 &\phi^N\psi^{NN} &\phi^B\psi^{BN} \\ 0 & 0 &\phi^3\alpha^3 &\dots &\phi^7\alpha^7 &0 &\phi^N\psi^{NB}+\iota &\phi^B\psi^{BB}+\iota \end{bmatrix}_{t} \begin{bmatrix} N_{1,t}|N_{1,t-1}\\N_{2,t}|N_{2,t-1}\\N_{3,t}|N_{3,t-1}\\N_{4,t}|N_{4,t-1}\\ \vdots\\ N_{7,t}|N_{7,t-1}\\ N_{7+,t}|N_{P,t-1}\\ N_{N,t}|N_{N,t-1}\\N_{B,t}|N_{B,t-1}\end{bmatrix}\)


where for year \(t\),

  • \(N_{j,t}\) is the number of pre-breeders for ages \(j\)=1,…,7

  • \(N_{7+,t}\) is the number of pre-breeders aged 7 or above never observed breeding

  • \(N_{N,t}\) is the number of mature non-breeders

  • \(N_{B,t}\) is the number of mature breeders

  • \(\phi^0_t\) is pup survival from birth to tagging at 1.5 months of age

  • \(\phi^1_t\) is apparent survival from tagging to the end of the first census year

  • \(\phi^j_t\) is apparent survival of seals aged \(j\)=2,…,7 between \(t\) and \(t\)+1

  • \(\phi^{7+}_t\) is apparent survival of pre-breeders never observed breeding between \(t\) and \(t\)+1

  • \(\phi^N_t\) is apparent survival of mature non-breeders between \(t\) and \(t\)+1

  • \(\phi^B_t\) is apparent survival of mature non-breeders between \(t\) and \(t\)+1

  • \(\alpha^j_t\) is the probability that a pre-breeder of age \(j\) starts breeding at that age between \(t\) and \(t\)+1

  • \(\psi^{B-}_t\) is the probabilities of being \(B\) in \(t\) and \(B\) or \(N\) in \(t\)+1, conditional on surviving

  • \(\psi^{N-}_t\) is the probabilities of being \(N\) in \(t\) and \(B\) or \(N\) in \(t\)+1, conditional on surviving

  • \(\varsigma_{F}\) is the proportion of females at birth

  • \(\iota_t\) is the number of breeding immigrants at \(t\)+1 per mature female at \(t\).


The IPM’s likelihood is \(\mathcal{L}(h,\phi,\alpha,\psi,p,y,n_{s}|N_{B},\iota)=\mathcal{L}_{MSCJS}(h|\phi,\alpha,\psi,p)\mathcal{L}_{sys}(N|\phi,\alpha,\psi,\iota)\mathcal{L}_{obs}(y|N_{B})\mathcal{L}_{c}(n_{s}|\beta,b,\lambda,\varphi,\zeta)\), and includes a multi-state Cormack-Jolly-Seber capture-mark-recapture model formulated as a Bayesian state-space model (Gimenez et al. 2007; Kéry and Schaub 2012b) with a state process and an observation process, denoted together as \(\mathcal{L}_{MSCJS}(h|\phi,\alpha,\psi,p)\). The model includes 12 states and represents survival and accession to breeding for pre-breeders (\(P\)) of ages 3 to 7, survival and transition of mature breeders (\(B\)) and for skippers or non-breeders (\(N\)), and two absorbing states, one for pre-bredeers older than 7 which are never seen breeding, and on for dead (\(D\)) or permanently emigrated females. The main data are the female encounter histories (\(h\)) as defined by categories specified in the observation process. The state process for marked individual \(i\), age \(j\), and year \(t\) is


\(z_{i,j,t+1} \sim \begin{pmatrix} 0 &\phi^1_{i,t} &0 &0 &0 &0 &0 &0 &0 &(1-\phi^1_{i,t}) \\ 0 &0 &\phi^2_{i,t} &0 &0 &0 &0 &0 &0 &(1-\phi^2_{i,t})\\ 0 &0 &0 &\phi^3_{i,t}(1-\alpha^3_{i,t}) &0 &0 &0 &\phi^3_{i,t}\alpha^3_{i,t} &0 &(1-\phi^3_{i,t})\\ 0 &0 &0 &0 &\ddots &0 &0 &\vdots &0 &\vdots \\ 0 &0 &0 &0 &0 &\phi^7_{i,t}(1-\alpha^7_{i,t}) &0 &\phi^7_{i,t}\alpha^7_{i,t} &0 &(1-\phi^7_{i,t}) \\ 0 &0 &0 &0 &0 &0 &\phi^7+_{i,t} &0 &0 &(1-\phi^{7+}_{i,t}) \\ 0 &0 &0 &0 &0 &0 &0 &\phi^B_{i,t}\psi^{BB}_{i,t} &\phi^B_{i,t}\psi^{BN}_{i,t} &(1-\phi^B_{i,t}) \\ 0 &0 &0 &0 &0 &0 &0 &\phi^N_{i,t}\psi^{NB}_{i,t} &\phi^N_{i,t}\psi^{NN}_{i,t} &(1-\phi^N_{i,t}) \\ 0 &0 &0 &0 &0 &0 &0 &0 &0 &1 \end{pmatrix}\)


where rows and columns are departure and arrival states between years \(t\) and \(t\)+1, including an absorbing state, \(P^{7+}\), for pre-breeders which are never observed breeding. The observation process \(Y_{i,t}\) is conditional on the true individual states, \(z_{i,t}\), in season \(t\) and comprises five types of female observations:

  • \(P^0\) - observed as newborn, subsequently tagged and released as a pup

  • \(P^{a}; a=1,...,7+\) - observed as pre-breeder of ages 1 to 7 or above (7+)

  • \(B\) - observed as breeder

  • \(N\) - observed as non-breeder

and an unobserved category, \(D\), which is dead or permanently emigrated. The observation process relates to the true states through the probabilities of detecting individuals observed as \(s\) at time \(t\), \(p^s_t\), and is formulated as

\(Y_{i,t+1}|Z_{i,t+1} \sim \begin{pmatrix} 0 &0 &0 &0 &1 \\ 0 &p^{P_1}_{t} &0 &0 &1-p^{P_1}_{t} \\ 0 &p^{P_1}_{t} &0 &0 &1-p^{P_1}_{t} \\ 0 &p^{P_1}_{t} &0 &0 &1-p^{P_1}_{t} \\ 0 &p^{P_2}_{t} &0 &0 &1-p^{P_2}_{t} \\ 0 &p^{P_2}_{t} &0 &0 &1-p^{P_2}_{t} \\ 0 &p^{P_2}_{t} &0 &0 &1-p^{P_2}_{t} \\ 0 &p^{P_2}_{t} &0 &0 &1-p^{P_2}_{t} \\ 0 &p^{P_2}_{t} &0 &0 &1-p^{P_2}_{t} \\ 0 &0 &p^{B}_{t} &0 &1-p^{B}_{t} \\ 0 &0 &0 &p^{N}_{t} &1-p^{N}_{t} \\ 0 &0 &0 &0 &1\end{pmatrix}\)


where rows are states \(P^a\), for \(a=1,..., 7+\), \(B\), \(N\), and \(D\), and columns are observation types in year \(t\). Note that detection probabilities are combined for ages 1 to 3 (\(p^{P_1}\)), which are less likely to be encountered, and for ages 4 to 7+ (\(p^{P_2}\)), to avoid parameter identifiability problems due to lack of data when fitting the models.


The matrix population model, with a system process and an observation process, \(\mathcal{L}_{sys}(N|\phi,\alpha,\psi,\iota)\mathcal{L}_{obs}(y|N_{BI})\), contributes to estimate population sizes and demographic rates, including immigration (\(\iota\)). Pup production data (\(y|N_{B}\)), which is observed with little error at SSB (Forcada et al. 2005), is equivalent to the observed breeding female population, and defines the observation process. For model fitting, I selected Poisson distributions to introduce demographic stochasticity, as follows.

\(N_{P_{1},t+1} \sim Poisson(\varsigma_{F}\phi^0_t\phi^1_tN_{BI,t})\)

\(N_{P_{2},t+1} \sim Poisson(\phi^2_tN_{2,t})\)

\(N_{P_{j},t+1} \sim Poisson(\phi^j_t(1-\alpha^j_t)N_{P_{j-1,t}}), \ {j=3,...,7+}\)

\(N_{R,t+1} \sim \sum_{j}Poisson(\phi^j_t\alpha^j_tN_{P_{j-1,t}}), \ {j=3,...,7}\)

\(N_{N,t+1} \sim \left[Poisson(N_{N,t}\phi^N_t\psi^{NN}_t)+Poisson(N_{BI,t}\phi^B_t\psi^{BN}_t)\right]\)

\(N_{B,t+1} \sim \left[Poisson(N_{N,t}\phi^N_t\psi^{NB}_t)+Poisson(N_{BI,t}\phi^B_t\psi^{BB}_t)+\sum_{j}Poisson(\phi_{j,t}\alpha^j_tN_{P_{j-1},t})\right]\)

\(N_{I,t+1} \sim Poisson\left[(N_{N,t}+N_{B,t})\iota_t\right]\)

Note that the observed pup production (\(y_t\)) results from the number of breeders, which includes new recruits from SSB, breeders already at SSB, and breeding immigrants from elsewhere. The distinction between \(N_B\) and \(N_{BI}\) allows for estimates of the fecundity (or pregnancy) rate, \(F_t\), which is the number of mature females (\(N_B+N_N\)) in year \(t\) which survive and breed in year \(t\)+1, excluding recruits and immigrants in year \(t\)+1, divided by the number of mature females surviving to year \(t\)+1; and importantly, estimates of number of immigrants, \(N_I\). This is achieved through the observation process as

\(y \sim Poisson(N_{BI,t})\)

\(N_{BI,t} = N_{B,t}+N_{I,t}\)

The last component of the likelihood, \(\mathcal{L}_{c}(n_{s}|\beta,b,\lambda,\varphi,\zeta)\), evaluates the proportion of the estimated SSB population, \(E([N_{t+1}|N_{t}])\), that represents a total daily female count at the SSB (denoted with superscript \(S\)). For this, daily counts obtained between November 1st and December 31st (\(n^S_k\) for day \(k\)) are modelled with a generalised additive mixed model (GAMM) with a Poisson-log-normal structure (Gelman and Hill 2006), which takes into account overdispersion. The more general model also considered autocorrelation of first order in female counts:

\(n_{k,t} \sim Poisson(\mu_{k,t})\)

\(log(\mu_{k,t}) = \beta_0+\sum_{l}b_l(d_k)\beta_l+\varphi\epsilon_{k,t-1}+\zeta_{k,t}\)

\(\zeta_{k,t} \sim N(0,\sigma^2_\varepsilon)\)

where \(d_k\) is the \(k\)th day of the survey season, \(\beta\) are unknown parameters (penalties), \(b\) is a vector of \(l\) spline basis functions, \(\varphi\) is an autocorrelation parameter with values between -1 and 1, \(\epsilon_{k,t-1}\) is the residual associated with observation \(k\) at day t-1, and \(\zeta_{k,t}\) an overdispersion parameter. The likelihood is \(\mathcal{L}_{c}(n_{s}|\beta,b,\lambda,\varphi,\zeta)\), where \(\lambda\) are smoothing parameters.

The combined likelihood produces derived quantities \(N^S_{M,t}⁄E(n^S_j)\) and \(N^S_{B,t}⁄E(n^S_j)\), where \(M=N+B\), \(N\) are mature non-breeding and \(B\) mature breeding females respectively; and \(E(n^S_j)\) is the predicted female count at SSB on the \(jth\) day. These multipliers of observed counts scale up estimates to total number of mature

\(\hat{N}_M=\sum_{c}n^c_j\hat{N}^S_M/E(n^S_j)\)

and breeding females

\(\hat{N}_B=\sum_{c}n^c_j\hat{N}^S_B/E(n^S_j)\)

at other locations, where \(n^c_j\) is the female count at colony \(c\) on day \(j\). Note that the second quantity \((\hat{N}_B)\) is equivalent to pup production, and the first quantity is the parameter of interest for the Red List assessment.

These equations assume that the dynamics of daily counts at SSB and other locations are highly correlated, which is accurate within Bird Island and across other locations, as previously shown (Forcada et al. 2023b). However, correlation between locations decays fast after the first six days following peak pupping, largely due to differences in haul out behaviour. This is shown in Fig.3 for female counts during season 2022.


**Figure 3**: Correlation of female counts at SSB and Maiviken (MVK; 105km E of SSB). Top left panel is the Pearson's correlation coefficient of predicted (from GAMM model) and observed counts at SSB during peak pupping, and at Maiviken from 0 to 20 days after peak, for 2014 to 2025. The top right panel shows in detail the maximum correlation of observed counts, with the regression of predicted counts at MVK on day 6 after peak against SSB at peak; counts at both locations were detrended to remove observed long-term trends. The bottom panels show the fit of Poisson GAM models with predictions rescaled to provide proportional predicted counts (with values between 0 and 1) for breeding females and territorial males during season 2022, when RPAS flights were conducted. Dots are observed counts, the black dashed vertical line is the estimated day of peak pupping, red and grey shaded areas correspond to negative and positive differences in counts between locations due to differences in seal haulout behaviour, and $\delta$ coefficients are predicted differences between SSB and MVK on days when RPAS counts were obtained.

Figure 3: Correlation of female counts at SSB and Maiviken (MVK; 105km E of SSB). Top left panel is the Pearson’s correlation coefficient of predicted (from GAMM model) and observed counts at SSB during peak pupping, and at Maiviken from 0 to 20 days after peak, for 2014 to 2025. The top right panel shows in detail the maximum correlation of observed counts, with the regression of predicted counts at MVK on day 6 after peak against SSB at peak; counts at both locations were detrended to remove observed long-term trends. The bottom panels show the fit of Poisson GAM models with predictions rescaled to provide proportional predicted counts (with values between 0 and 1) for breeding females and territorial males during season 2022, when RPAS flights were conducted. Dots are observed counts, the black dashed vertical line is the estimated day of peak pupping, red and grey shaded areas correspond to negative and positive differences in counts between locations due to differences in seal haulout behaviour, and \(\delta\) coefficients are predicted differences between SSB and MVK on days when RPAS counts were obtained.

As seasonal count data were also available for Maiviken (MVK), I incorporated differences in haul out patterns between SSB and MVK to make observed counts at or near MVK scalable and comparable according to the ratios \(N^S_{M,t}⁄E(n^S_j)\) or \(N^S_{B,t}⁄E(n^S_j)\), which relate specifically to SSB counts. For this, I obtained predicted counts at both locations from Poisson GAMs fited to daily data, and predictions were rescaled to range between 0 and 1; \(z_i=(n_i-min(n_i))/(max(n_i)-min(n_i))\).

To obtain a corrected MVK count on day \(i\), I assume a difference of size \(\kappa_i\) in \(E(n_i)\) between locations. At the transformed scale,

\(\kappa=(n^M_i\delta_i)/z^M_i\), where \(\delta^-_i\) is \(z^M_i-z^S_i\), when \(z^M_i>z^S_i\), and \(\delta^+_i\) is \(z^S_i-z^M_i\) when \(z^S_i>z^M_i\). See for example lower panels of Fig. 3. Thus, a rescaled MVK count with respect to SSB, is \(n^M_iK_i\), where

\(K_i=1\), if \(z^M_i=z^S_i\)

\(K_i=1+\frac{\delta^-_i}{z^M_i}\), if \(z^M_i>z^S_i\)

\(K_i=1-\frac{\delta^+_i}{z^M_i}\), if \(z^M_i<z^S_i\).

This adds prediction uncertainty to counts away from SSB, more likely to follow a MVK haul out pattern, and for counts closer to SSB, some of the uncertainty can be reduced by obtaining and average corrected count at colony \(c\), \(\overline{n}^c_j\), weighted by the inverse of the distance in km from \(c\) to SSB and to MVK. This is obtained as

\(\overline{n}^c_j=\frac{(w_{1,j}n^c_j+w_{2,j}n^c_jK_j)}{\sum{w_j}}\),

where \(w_{1,j}=1/d^{SSB}_j\) and \(w_{2,j}=1/d^{MVK}_j\), with \(d\) as distance in km. Thus, a count closer to SSB will assume a more similar daily haul out pattern to this location, with \(\overline{n}^c_j\rightarrow n^c_j\), whereas a count closer to MVK will have a stronger influence of \(K_j\), with \(\overline{n}^c_j\rightarrow n^c_jK_j\).

It follows that an estimate of total number of mature females at a survey area is
\(\hat{N}_{F_M}=\frac{\hat{N}^S_M}{E(n^S_j)}\sum_{c}\frac{w_1n^c_j+w_2n^c_jK_j}{w_1+w_2}\),

and to account for pre-breeders and obtain an estimate of total number of females at a survey area I use
\(\hat{N}_{F_{All}}=\left(1+\frac{\sum_{k}\hat{N}^S_{P_{k}}}{\hat{N}^S_M}\right)\frac{\hat{N}^S_M}{E(n^S_j)}\sum_{c}\frac{w_1n^c_j+w_2n^c_jK_j}{w_1+w_2}\),

where \(\hat{N}^S_{P_{k}}\) is the estimated number of pre-breeders of age \(k\), for \(k=1,...,7+\), at SSB.


2.3.1.1.1.2 Males

The male IPM has a combined age-stage matrix model with pre-territorial males (\(P\)) of ages 1 to 11+ (ages 11 or older), and mature states including Territorial (\(T\)) and mature non-territorial (\(N\)). The lowest observed age at first territory tenure is 7, and first-time territorials older than 11 are seldom observed. The expected number of individuals in year \(t\)+1 is


\(E([N_{t+1}|N_{t}])=E\left( \begin{bmatrix} N_{1,t+1}|N_{1,t}\\N_{2,t+1}|N_{2,t}\\ \vdots\\ N_{7,t+1}|N_{7,t}\\ \vdots\\ N_{11,t+1}|N_{11,t}\\ N_{N,t+1}|N_{N,t}\\N_{T,t+1}|N_{T,t}\end{bmatrix}\right)=\)


\(\begin{bmatrix} \varsigma_{M}\phi^0\phi^1 &0 &0 &0 &0 &0 &0 &0\\0 &\phi^2 &0 &0 &0 &0 &0 &0\\0 &0 &\ddots &0 &0 &0 &0 &0\\ 0 &0 &0 &\phi^7(1-\alpha^7) &0 &0 &0 &0\\ 0 &0 &0 &0 &\ddots &0 &0 &0\\0 &0 &0 &0 &0 &\phi^{11}(1-\alpha^{11}) &0 &0\\ 0 &0 &0 &0 &0 &0 &\phi^N\psi^{NN} &\phi^T\psi^{TN}\\0 &0 &0 &\phi^7\alpha^7 &\dots &\phi^{11}\alpha^{11} &\phi^N\psi^{NT} &\phi^T\psi^{TT}\end{bmatrix}_{t}\begin{bmatrix} N_{1,t}|N_{1,t-1}\\N_{2,t}|N_{2,t-1}\\\vdots\\ N_{7,t}|N_{7,t-1}\\ \vdots\\N_{11,t}|N_{11,t-1}\\ N_{N,t}|N_{N,t-1}\\N_{T,t}|N_{T,t-1}\end{bmatrix}\)


where for year \(t\),

  • \(N_{j,t}\) is the number of pre-territorial males for ages \(j\)=1,…,11

  • \(N_{N,t}\) is the number of territorial males skipping a breeding season

  • \(N_{T,t}\) is the number of territorial males

  • \(\phi^0_t\) is pup survival from birth to tagging at 1.5 months of age

  • \(\phi^1_t\) is apparent survival from tagging to the end of the first census year

  • \(\phi^j_t\) is apparent survival of seals aged \(j\)=2,…,11, between \(t\) and \(t\)+1

  • \(\phi^N_t\) is apparent survival of territorial males skipping a breeding season between \(t\) and \(t\)+1

  • \(\phi^T_t\) is apparent survival of territorial males between \(t\) and \(t\)+1

  • \(\alpha^j_t\) is the probability that a male of age \(j\) starts holding territory at that age between \(t\) and \(t\)+1

  • \(\psi^{T-}_t\) is the probabilities of being \(T\) in \(t\) and \(T\) or \(N\) in \(t\)+1, conditional on surviving

  • \(\psi^{N-}_t\) is the probabilities of being \(N\) in \(t\) and \(T\) or \(N\) in \(t\)+1, conditional on surviving

  • \(\varsigma_{M}\) is the proportion of females at birth

  • \(\iota_t\) is the the proportion of immigrants among territorial males at \(t\)+1.


The male IPM’s likelihood is \(\mathcal{L}(h,\gamma,\phi,\alpha,\psi,p,y,n_{s}|N_{B},\iota)=\mathcal{L}_{MSCJS}(h|\phi,\alpha,p)\mathcal{L}_{MSJS}(h|\phi,\gamma,\psi,p)\mathcal{L}_{sys}(N|\phi,\alpha,\psi,\iota)\mathcal{L}_{obs}(y|N_{T})\mathcal{L}_{c}(n_{s}|\beta,b,\lambda,\varphi,\zeta)\). It is similar to the female IPM likelihood, but the capture-mark-recapture component is split into separate likelihoods given the different data types (PIT-tag vs genetic recaptures), sample size, and short temporal overlap of two data sets. Pups with implanted PIT-tags at birth represent largely pre-territorial component, and are used to model survival and accession to territories, \(\phi^j_{j,t}\) and \(\alpha^j_t\), with likelihood \(\mathcal{L}_{MSCJS}(h|\phi,\alpha,p)\); and territorial males captured and recaptured through genetic sampling contribute to estimate their survival (\(\phi^{T,N}_{t}\)), territory tenure (\(psi^{-T}_t\)) and temporary emigration rates (\(psi^{-N}_t\)), with likelihood \(\mathcal{L}_{MSJS}(h|\phi,\gamma,\psi,p)\).

In the MSCJS model the state process equation for marked pre-territorial individual \(i\) of age \(j\) in year \(t\), is formulated as


\(z_{i,t+1} \sim \begin{pmatrix} 0 &\phi^1_{i,t} &0 &0 &0 &0 &0 &0 &0 &(1-\phi^1_{i,t}) \\ 0 &0 &\phi^2_{i,t} &0 &0 &0 &0 &0 &0 &(1-\phi^2_{i,t})\\ 0 &0 &0 &\ddots &0 &0 &0 &0 &0 &\vdots\\0 &0 &0 &0 &\phi^7_{i,t}(1-\alpha^7_{i,t}) &0 &0 &0 &\phi^7_{i,t}\alpha^7_{i,t} &(1-\phi^7_{i,t}) \\ 0 &0 &0 &0 &0 &\ddots &0 &0 &\vdots &\vdots \\ 0 &0 &0 &0 &0 &0 &(1-\phi^{10}_{i,t}\alpha^{10}_{i,t}) &0 &\phi^{10}_{i,t}\alpha^{10}_{i,t} &(1-\phi^{10})\\ 0 &0 &0 &0 &0 &0 &0 &\phi^{11}_{i,t} &0 &(1-\phi^{11}_{i,t}) \\ 0 &0 &0 &0 &0 &0 &0 &0 &\phi^T_{i,t} &(1-\phi^T_{i,t}) \\ 0 &0 &0 &0 &0 &0 &0 &0 &0 &1 \end{pmatrix}\)


where rows and columns are departure and arrival states between years \(t\) and \(t\)+1, assuming that at age 11 all males which survive are territorial. The observation process \(Y_{i,t}\) conditional on true individual states (\(z_{i,t}\)) comprises three types of observations:

  • \(P^0\) - observed as newborn, subsequently tagged and released as a pup

  • \(P^{a}; a=1,...,7+\) - observed as pre-territorial of ages 1 to 11

  • \(T\) - observed as territorial.

and an unobserved category, \(D\), which is dead or permanently emigrated. The observation process includes the probabilities of detecting individuals observed as \(s\) at time \(t\), \(p^s_t\), and is formulated as


\(Y_{i,t+1}|Z_{i,t+1} \sim \begin{pmatrix} 0 &0 &0 &1\\ 0 &p^{P_1}_{t} &0 &1-p^{P_1}_{t}\\ 0 &\vdots &0 &\vdots\\0 &p^{P_{11}}_{t} &0 &1-p^{P_{11}}_{t}\\ 0 &0 &p^{T}_{t} &1-p^{T}_{t}\\0 &0 &0 &1\end{pmatrix}\)

where rows are states \(P^0\) (pup), \(P^a\), for \(a=1,..., 11\), \(T\), and \(D\), and columns are observation types in year \(t\).

As the elapsed time between birth and accession to territory can be 10 years or more, observations of pre-territorials becoming territorials over 25 years of study were only a small fraction of the total number of PIT-tagged male pups. This limited parameter estimability, statistical power and model selection, making IPM fitting inefficient. As a valid alternative (Besbeas et al. 2002), I used informative priors for \(ϕ^j_t\) and \(\alpha^j_t\) in the IPMs, obtained from a simplified MSCJS recruitment model (Pradel and Lebreton 1999) fitted independently to the PIT-tag data set. For this, I combined detection probabilities for pre-territorial of ages 1 to 11 (\(p^{P_j}=p; j=1,...,11\)), which are mostly unobservable, and set \(\alpha^j, j=7,...,11\), constant for the \(j\)th age.

The MSJS likelihood, \(\mathcal{L}_{MSJS}(h|\phi,\gamma,\psi,p)\), is based on a multi-state Jolly-Seber model with hierarchical state-space formulation and parameter-expanded data augmentation (Royle and Dorazio 2009; Kéry and Schaub 2012a). The model considers only territorial males. The total number of individuals alive during the study is \(N^s\), and \(b_t\) is the probability that a new territorial male enters the population in year \(t\), so that the number of new territorial males in a year is

\(B_t=N^sb_t\).

For data augmentation (Kéry and Schaub 2012a), I consider \(M\) individuals, of which \(N^s\) is a subset, and \(b_t\) is derived as a restricted dynamic occupancy model (Royle, Dorazio, and Link 2007). This relies on \(\gamma_t\), which is the probability that a male from \(M\) holds territory for the first time in year \(t\), so that

\(E(B_t)=\prod_{i=1}^{t-1}(1-\gamma_i)\gamma_t\),

and

\(N^s=\sum_{t}B\).

The state process equation for male \(i\) and year \(t\), and state \(r\) is formulated as


\(z_{i,t+1} \sim \begin{pmatrix} 1-\gamma_t &\gamma_t &0 &0\\ 0 &\phi^T_{i,t}\psi^{TT}_{i,t} &\phi^T_{i,t}\psi^{TN}_{i,t} &1-\phi^T_{i,t}\\ 0 &\phi^N_{i,t}\psi^{NT}_{i,t} &\phi^N_{i,t}\psi^{NN}_{i,t} &1-\phi^{N}_{i,t}\\ 0 &0 &0 &1\end{pmatrix}\)

where the states in rows and columns are \(P\) (not yet territorial in the population), \(T\), \(N\), and \(D\) (Dead or permanently emigrated), between years \(t\) and \(t\)+1. The observation process is


\(Y_{i,t+1}|Z_{i,t+1} \sim \begin{pmatrix} 0 &1\\ p^T_t &1-p^T_t\\ 0 &1\\ 0 &1\end{pmatrix}\)

where \(p^T_t\) is the recapture probability of territorial males and non-territorial males are unobserved, as in a temporary emigration model (Schaub et al. 2004).

The number of new territorial males each year (\(B_t\)) confounds males born at SSB starting to hold territory with immigrants that never held a territory at SSB before, and the two components are separated with the demographic part of the IPM (below), which combines both male capture-recapture data sets.

The matrix population model, \(\mathcal{L}_{sys}(N|\phi,\alpha,\psi,\iota)\mathcal{L}_{obs}(y|N_{T})\), with data on male pup production, \(\varsigma_{M}N_{0,t}\), which defines the observation process, contributes to estimate population sizes and demographic rates, including immigration (\(\iota\)). For model fitting, I select Poisson distributions which introduce demographic stochasticity. For the pre-territorial component, the equations are

\(N_{P_{1,t+1}} \sim Poisson(\varsigma_{M}\phi^0_t\phi^1_tN_{0,t})\)

\(N_{P_{j,t+1}} \sim Poisson(\phi^j_t N_{P_{j,t+1}}) , \ {j=2,...,6}\)

\(N_{P_{k,t+1}} \sim Poisson(\phi^k_t(1-\alpha^k_t)N_{P_{k-1,t}}) , \ {k=7,...,11}\)

For the territorial data set, the MSJS model directly produces estimates of \(N_{N,t+1}\), \(N_{T,t+1}\), and \(B_t\), which allows for an estimate of immigrants as

\(N_{I,t+1} \sim B_t - \sum_{j}Poisson(\phi^j_t\alpha^j_tN_{p_{j-1,t}})\)

The fraction of the estimated representative study population, \(E\left([N_{t+1}|N_t]\right)\), that is a territorial male count (\(n^S_j\)) obtained between November 1st and December 31st, is modelled with a Poisson GAMM, exactly as in the female IPM. It is the last component of the male IPM likelihood, \(\mathcal{L}_{c}(n_{s}|\beta,b,\lambda,\varphi,\zeta)\) and results in derived quantities \(N^S_{M,t}⁄E(n^S_j)\) and \(N^S_{T,t}⁄E(n^S_j)\), where \(M\) and \(T\) are mature (\(T+N\)) and territorial (\(T\)) males respectively; and \(E(n^S_j)\) which is the predicted male count at SSB on the \(j\)th day.

These are Used as multipliers of observed counts at other locations to obtain estimates of total number of mature males

\(\hat{N}_M=\sum_{c}n^c_j\hat{N}^S_M/E(n^S_j)\)

and territorial males

\(\hat{N}_T=\sum_{c}n^c_j\hat{N}^S_T/E(n^S_j)\),

where \(n^c_j\) is the male count at colony \(c\) on day \(j\).

The estimate of mature males at a study area adjusting for differences in haul out between SSB and MVK is

\(\hat{N}_{M_M}=\frac{\hat{N}^S_M}{E(n^S_j)}\sum_{c}\frac{w_1n^c_j+w_2n^c_jK_j}{w_1+w_2}\),

and an absolute estimate of total number of males at a study area is obtained adjusting for pre-territorial males,

\(\hat{N}_{M_{All}}=\left(1+\frac{\sum_{k}\hat{N}^S_{P_{k}}}{\hat{N}^S_M}\right)\hat{N}_{M_M}\),

where \(\hat{N}^S_{P_{k}}\) is the estimated number of pre-territorials of age \(k\), for \(k=1,...,11\), at SSB. (See female section for the derivation of \(\overline{n}^c_j\).

2.3.1.1.2 IPMs implementation, assessment, and selection

I closely follow methods of Forcada et al. (2023b) and make the assumption of independence between data sets. As in many IPMs, this is not met when for example, I use the encounter histories \(h\) of individuals which are also counted. However, based on recent experimental evidence (Abadi et al. 2010; McDonald et al. 2016; Weegman et al. 2021), this is unlikely to be a source of bias and affect model estimation and performance, provided that model assumptions are met (see Riecke et al. 2019). These assumptions have already been evaluated for earlier versions of the the same data sets in model components of our IPMs which have been used independently elsewhere (Forcada, Trathan, and Murphy 2008; Forcada and Hoffman 2014). For example, the effects of tag-loss on survival. However, when fitting the models, I used predictive checks to assessed the ability of the models to generate data that were consistent with the observed data (Gelman et al. 2002) (data not presented), and in particular in relation to variation between estimates of survival and emigration derived from the IPM.

I use Markov Chain Monte Carlo (MCMC) methods in BUGS language and program JAGS (Plummer 2003), run from \(\texttt{R}\) (v4.5.0 R_Core_Team 2024) using packages \(\texttt{jagsUI}\) (Kellner 2024) and \(\texttt{rjags}\) (Plummer 2024) to fit the IPMs. To specify priors, I use Normal distributions with numbers of pre-breeders, breeders and mature non-breeders observed in the first year for the initial population sizes. For the capture-recapture part, model fit was initially assessed with frequentist multistate goodness-of-fit testing (Pradel, Wintrebert, and Gimenez 2003) using \(\texttt{R2ucare}\) (Gimenez et al. 2018). Selection of best models was subsequently based on WAIC (Gelman, Hwang, and Vehtari 2014; Watanabe 2010), and the IPMs were fitted with the best mark-recapture modelling options.

I selected temporal variation in parameters as much as possible to reflect the population state at the time of data collection and obtain accurate abundance estimates; but it was not always possible when data were sparse as in the male mark-recapture. For the female IPM, I used independent temporal random effects for \(\phi_j\) and \(\alpha_j\) for ages 1 to 7 and 3 to 7 respectively, \(\psi^{NB}\), \(log(\iota)\), \(p^{P1}\), \(p^{P2}\), \(p^B\) and \(p^N\) (e.g. \(logit(\psi^{NB}_{i,t})=\mu^{^NB}_i+\varepsilon^{NB}_{i,t}\) , where \(\varepsilon^{NB} \sim Normal(0,\sigma^2_{NB})\), with diffuse Uniform priors for \(\sigma^2_*\). I used correlated temporal random effects for \(\phi_B\), \(\phi_N\) and, \(\psi_{BN}\), where \(\varepsilon^{\theta}_{i,t} \sim MVNormal(0,\Sigma_{\theta,t})\) and \(\Sigma_{\theta,t}\) is a temporal variance-covariance matrix describing the temporal variance of each parameter (\(\theta\)) and the temporal covariance among parameters (Link and Barker 2005). For \(\Sigma_{\theta,t}\) I used priors from an inverse Wishart distribution (Gelman and Hill 2006) with scale matrix \(\Omega\) as identity matrix and for degrees of freedom. For the genetic mark-recapture part of the male IPM, I chose the best combination of \(\gamma\), \(\phi\), and \(\psi\), either with temporal random effects, correlated temporal random effects, or constant values over time based on the lowest WAIC of a fitted set of models with different parameter combinations (data not presented). For juveniles and pre-territorials, \(\phi_j\) and \(\alpha_j\) were selected as informative Beta priors with shape parameters obtained from an independent best fit of a multistate mark-recapture model of the PIT-tag data set. The priors incorporated temporal variation as observed in the data record, reflecting uncertainty due to environmental stochasticity at the time of the survey counts.

The GAMMs of daily counts in the IPMs were implemented using methods of Wood (2016). Thus, I used multivariate Normal priors on the smooth coefficients in the JAGS models, and the smoothing parameters and smoothing penalty matrices directly specified the prior multivariate normal precision matrix as input data. I implemented preliminary model selection based on WAIC to test for overdispersion and autocorrelation in daily counts before incorporating the best Poisson models in the IPMs.

In all IPM analyses, I used 200,000 iterations of three Markov chains and discarded the first 50,000 samples of each chain as burn-in phase, thinning the remainder to every 10th sample. I assessed chain convergence visually using trace plots, through the mixing of the chains and sample autocorrelation plots, and using the \(\hat{R}\) potential scale reduction factor statistic of less than 1.05 to retain posterior distribution samples (Gelman and Shirley 2011; Vehtari et al. 2021).

2.3.1.2 South Orkney Islands

2.3.1.2.1 Population dynamics

I evaluated population trends at Signy Island using a log-linear model with a discrete-time stochastic Gompertz formulation (Dennis and Taper 1994; Lebreton and Gimenez 2013) of the observed fur seal count \(y_t\) for year \(t\), in response to a year-lagged count, \(y_{t-1}\), with coefficient \(ν\). To accommodate positive skewness from numerous low counts, I used a negative binomial error distribution (e.g., Davis and Wu 2009), with \(y_t \sim NB(r,p_t)\), where \(r\) and \(p_t\) were the dispersion and success rate parameters respectively, so that \(p_t=r/(r+E[y_t])\).

Expressed in natural log scale, \(x_t=log E[y_t]\), the Gompertz Negative Binomial Model (GNBM) was formulated as

\(log (\frac{x_t}{x_{t-1}})=\upsilon+(1+\nu)y^T_{t-1}+\varepsilon\)

where intercept \(\upsilon\) measures intrinsic growth rate, equivalent to \(log(\lambda)\) when \(E[y]=1\), \(\nu\) is a density-dependence coefficient, and \(\varepsilon\) and error term.

2.3.1.2.2 Synchrony between South Orkneys and South Georgia

To formally test for synchrony between locations, I use a mixed-effects formulation (GNBMM) with two independent random effects (Grosbois et al. 2009; Lahoz-Monfort et al. 2011; Santin-Janin et al. 2014). A time-specific random effect \(\delta_t \sim N(0,\sigma^2_{\delta})\) which models common synchronous variation across sites; and a separate time-specific random effect for each site \(s\), \(\epsilon_{t,s} \sim N(0,\sigma^2_\epsilon)\) which accounts for additional and asynchronous variation. The effects are incorporated as

\(log (\frac{x_{s,t}}{x_{s,t-1}})=\upsilon_s+(1+\nu_{s,0})y^T_{s,t-1}+\delta_t+\epsilon_{t,s}\),

Here, \(\sigma^2_\delta+\sigma^2_{\epsilon_s}\) is between-year variance unexplained by density-dependence. The variances are equally partitioned between sites in the shared term (\(\sigma^2_\delta\)), but not between unshared components (\(\sigma^2_{\epsilon_s}\)). The fraction of between-year variance accounted by the shared component is the intra-class correlation,

\(ICC_s=\sigma^2_\delta/(\sigma^2_{\delta}+\sigma^2_{\epsilon_s})\),

which measures the synchrony of SI with BI, and BI with SI. When the shared (synchrony) component is large relative to the unshared component, \(ICC_s\) approaches 1, and tends to zero otherwise.

2.3.1.2.3 Connection between Signy Island and South Georgia

To further examine the association of the population abundance at Signy Island (SI) with the abundance at Bird Island (BI), I used a cross-correlation analysis function \(R_{y_{BI},y_{SI}}(τ)=E\left[y_{BI,t-τ},\overline{y}_{SI,t}\right]\), where \(τ\) is a difference in years, to uncover lagged correlations between both time-series. Based on these, I then investigated the population effects of BI, proxy for South Georgia, on temporary immigration at SI at appropriate time-lags (\(\tau\)). For that, I augmented the basic GNBM as

\(log (\frac{x_t}{x_{t-1}})=\upsilon+(1+\nu_0)y^T_{SI,t-1}+\nu_1y^T_{BI,t-\tau}+\varepsilon\),

where \(\nu_0\) is a density-dependence coefficient, and \(\nu_1\) is a fixed-effect of the BI-SSB mature female abundance at time-lag \(\tau\).

2.3.1.2.4 Model fitting, selection, and assessment

I fitted Gompertz negative binomial models using a state-spece framework and with the same software and approach as the IPMs described in the previous section. To specify priors, I selected \(r_s \sim U(0,50)\), where \(U\) is for Uniform distribution, as dispersion parameter for each site \(s\) (SI or BI). In single site GNBMs, for intrinsic growth rate I used \(\upsilon \sim N(10,10^{-2})\), where \(N\) is for Normal distribution; and in two-site GNBMs, for site \(s\) I used \(\upsilon_s \sim N[\bar\upsilon,U(0,2)^2]\) and \(\bar\upsilon \sim N(10,10^{-2})\). In single-site GNBMs, priors for density-dependence effects were \(\nu \sim N[0,U(0,2)^2)]\); and for density-dependence and lagged (\(\tau\)) BI effects on SI, these were \(\nu_1 \sim N[0,U(0,2)^2)]\) and \(\nu_2 \sim N(0,10^{-2})\). In single-site GNBMs, system process variation priors were \(\sigma^2_\varepsilon \sim \Gamma^{-1}(1^{-1},1^{-2})\) (where \(\Gamma\) is for Gamma distribution); and for two-site, synchrony GNBMMs, \(\sigma^2_{\varepsilon_s} \sim \Gamma^{-1}(1^{-1},1^{-1})\) for each \(s\), and \(\sigma^2_{\delta} \sim \Gamma^{-1}(1^{-1},1^{-1})\).

I used 350,000 iterations of four Markov chains and discarded the first 250,000 samples of each chain during the burn-in phase, thinning the remainder to every 25th sample, and thus producing 16,000 posterior simulations for inference. I assessed chain convergence visually using trace plots, through the mixing of the chains and sample autocorrelation plots, and using the \(\hat{R}\) potential scale reduction factor statistic of less than 1.05 to retain posterior distribution samples (Gelman and Shirley 2011; Vehtari et al. 2021). I assessed the ability of the models to generate data that were consistent with the observed data using posterior predictive checks (Gelman et al. 2002).

2.3.1.2.5 Scales and patterns of fluctuation in population time-series

To discriminate fluctuations from long-term trends, I investigated scales and patterns in population fluctuations across locations using a wavelet decomposition of the BI-SSB mature female abundance and SI population index. This highlighted short-lived (transient) dynamics, which combined led to mean synchrony across locations. I used wavelet power analyses \(W_x(f,\tau)\) for estimates at each location (\(x\)), and a cross-wavelet of both locations \(W_{x_1,x_2}(f,\tau)\), using R package \(\texttt{biwavelet}\) [v0.20.22; Gouhier, Grinsted, and V. (2024)]. The cross-wavelet helped identifying cyclical patterns in both time-series, how these varied across different time scales, and revealing common patterns of gradual change (Cazelles et al. 2008), and how these patterns differed from long-term declines.

2.3.2 IUCN Red List criteria selection and assessment

For assessments purposes, currently distinct genetic groups are considered subpopulations as defined by the IUCN; i.e., “geographically or otherwise distinct groups in the population between which there is little demographic or genetic exchange (typically one successful migrant individual or gamete per year or less)” (IUCN 2024). Here, I consider the four identified subpopulations under this definition, with the caveat that the genetic identity of the South Sandwich Islands remains to be ascertained.

Major observed subpopulation declines are environmentally driven and are ongoing, but the specific causes remain uncertain and exacerbated by other effects. Becasue the declines are likely to be irreversible in the near future, largely due to poor mitigation of climate change and other anthropogenic effects, my main focus was to evaluate status under IUCN Red List Criteria A1-A4, with special emphasis on criterion A2. This is defined as “population reduction observed, estimated, inferred, or suspected in the past where the causes of reduction may not have ceased or may not be understood or may not be reversible” (IUCN 2024). Other criteria were less likely, due to best evidence including the large size of the global population of the species, and absence of perceived significant change in range size or area of occupancy.

The assessment thus focused on population size reduction to assign Red List category. It was based on the mature population only, over the longer of 10 years or three generations. Due to major differences in life history between sexes, the mature population has a much higher proportion of females than males. According to IUCN guidelines, mature individuals are capable of reproduction, but “those who will never produce new recruits should not be counted” (IUCN 2024). Pre-territorial males, while physiologically able to reproduce from ages 3 to 4 (Payne 1979), they only start holding breeding territories when they are at least 7, at near physical maturation, but the mean age at first territory is between 9 and 11, depending on densities of individuals (Forcada 2021).

2.3.2.1 Generation length

This metric was evaluated as the average age of parents in the population (Coale 1972; IUCN 2024),

\(G=\frac{\sum_{x}x l_x m_x)}{\sum_{x}l_x m_x)}\),

where \(x\) is age and summations are from age 0 to the last observed age of reproduction, which is 24 for females and 18 for males; \(m_x\) is the fecundity at age; and \(l_x\) is survivorship up to age \(x\).

For this analysis, I used the posterior simulations of vital rates from the fitted IPMs to generate simulations of their core matrix population models. I then used age-from-stage decomposition methods (Caswell 2001) as implemented in \(\texttt{R}\) package \(\texttt{Rage}\) (O. R. Jones et al. 2022) to generate life table simulations for each sex, and derived generation length estimates.

As there is significant generation length variation in response to the environment (e.g., Forcada et al. 2023b), I used pre-disturbance vital rate values, corresponding to the period immediately prior to the population collapse of Bird Island in 2009. Additionally, because generation length varies between sexes, \(G\) was weighted by the relative frequency (abundance) of mature seals of each sex, obtained directly from the fitted IPMs.

2.3.2.2 Population change and reduction

Based on the observed population trajectories at Bird Island (Forcada et al. 2023b) and at the study areas of the South Shetland Islands (e.g., Figure3 in Krause et al. (2022)), I assumed an exponential decline model to calculate reduction of the subpopulations in the South Atlantic region, and of the total population. For this, I considered two populations, \(N_1\) which was population size in the first year (\(t_1\)), and \(N_T\) which was the size in the last year (\(t_T\)), and the observed period \(T-1\) was the number of years between the first and last observation, \(t_T-t_1\). The observed change, \(N_T/N_1\) was expressed as annual change as

\(AC=\left(\frac{N_T}{N_1}\right)^{\frac{1}{T-1}}\),

or as a percent annual change, as \(\%AC=100(1-AC)\). Population reduction was evaluated for an interval equivalent to three times the generation length (\(3G\)), and was estimated as

\(R=1-\left(\frac{N_T}{N_1}\right)^{\frac{3G}{T-1}}\).

I based calculations of population reduction under IUCN Red List Criteria A1-A4 on the Criterion A tool (https://www.iucnredlist.org/resources/criterion-a). However, I did it separately to include uncertainty in population estimates, generation times, and pre-assessment annual change rates (see below), which is not possible with the tool.

Because a \(3G\) length interval could precede and/or outlast the interval between the first and last population estimates available for the assessment, I selected the years of data from each breeding location which more closely matched the assessment period. When this was not possible but there were two separate population estimates, I obtained the population three generations ago as

\(N_{3G}=N_1/AC^{t}\),

where \(AC\) was the annual change over the number of years (\(t\)) between the beginning of the assessment period and the year when \(N_1\) was obtained. This \(AC\) was not always the same as the estimated \(AC\) between \(t_1\) and \(t_T\) (above), and it was therefore obtained from the best additional data available for each breeding location. For example, at South Georgia it was obtained from the geometric mean estimates of \(log\lambda\) derived from the IPMs. For other locations in Table 1, pup production values preceding, or soon after, the start of the assessment period were available, and this was less of an issue.

The current population (\(N_C\)), if the assessment year was later than \(t_T\), was obtained as

\(N_C=N_TAC^{t}\),

where \(AC\) was evaluated over the period between \(T\) and the assessment year (\(AY\)).

For this analysis, the assessment year was 2025, which gave an assessment period of nearly 26 years, starting in 1999, and was equivalent to three times the estimated generation length of 8.6 years (see Results section, below). For South Georgia, reduction was evaluated with data derived from this analysis and the published record. For other breeding locations, I used pup production/count data available in the published and unpublished records (Table 1).

2.3.2.2.1 Adjustment of pup counts to mature, and total population size

Only pup counts, representative of total pup production, were available from one or more breeding locations of South Shetland Islands, Bouvetøya, and East Antarctica (Table 1). Estimates of mature individuals were thus obtained from local pup production values (\(n^P\)) and the ratio of pup production and mature individual abundance estimates from the fitted IPMs. This was obtained as

\(\hat{N^M}=n_P\hat{\pi_M}\)

where the multiplier \(\hat{\pi_M}\) was

\(\hat{\pi_M}=\sum_{t=1}^{T}\left(\frac{N_t^{BI}+N_t^{NB}+N_t^T+N_t^{NT}}{N_t^{BI}}\right)T^{-1}\)

where in year \(t\), \(N^{BI}\) are breeding females including immigrants (see the female IPM formulation), \(N^{NB}\) are mature females skipping breeding, \(N^{T}\) are territorial males, and \(N^{NT}\) are mature non-territorial males. Averaged over a number of years, the estimate of \(\hat{\pi}\) reflects temporal variation, which was incorporated in extrapolations of pup counts. To incorporate further uncertainty, I used the IPM simulations of these parameters together with log-normal deviates of pup production based on their mean and variance, \(\hat{Var}(n_P)\), when the variance was available. Otherwise, I only used the mean (\(\overline{n}_P\)).

An adjustment from pup counts to total population was obtained as

\(\hat{\pi_T}=\sum_{t=1}^{T}\left(\frac{N_t^{BI}+N_t^{NB}+N_t^{PB}+N_t^T+N_t^{NT}+N_t^{PT}}{N_t^{BI}}\right)T^{-1}\)

where \(N^{PB}\) was the abundance of female pre-breeders and \(N^{PT}\) the abundance of pre-territorial males.

Previously, extrapolations from pup counts to total population (e.g., Hofmeyr et al. 2005; Krause et al. 2022) were based on adjustments derived by Payne (1979) from unrepresentative and incomplete data, and grossly overestimated total population size (Forcada et al. 2023b), especially of males. Additionally, these multipliers adjusted for total population, as opposed to mature population only, which was the objective of the assessment. Here, such old adjustments are deemed unreliable, as discussed below, and I derived robust adjustments from the IPMs, relevant to the assessment period.

2.3.2.2.2 South Georgia population projection

The total mature population size of South Georgia (\(N_1^M\)) was obtained from two comprehensive helicopter surveys of all the breeding locations conducted in 2007/09 (\(t_1\)), and the data are available at (Forcada et al. 2023a). The population size in 2022 (\(t_T\)) was only available for a few representative locations (Fig. 2a) and, therefore, I derived a population total using a projection model. The model assumed that abundance of each colony \(c\) surveyed in \(t_1\) changed at an annual rate \(\overline{AC}_c\), which was an average of the \(AC_s\) of the five colonies surveyed in 2007/9 and 2022, weighted by the distance (\(d\)) from \(c\) to each colony \(s\);

\(\overline{AC}_c=\frac{\sum_{s}d_{cs}AC_s}{\sum_{s}d_{cs}}\).

This ensured that regional differences in annual change where taken into account, and a colony \(c\) would change at an annual rate which was more similar to a rate observed at the \(s\) locations which were geographically closer.

Total mature abundance was obtained as

\(\hat{N}_T^M=\sum_{c}\hat{N}_{1,c}^M\hat{AC}_c^{(T-1)}\).

where the uncertainty of the estimated \(\overline{AC}_c\) reflected the uncertainty of the IPM simulations used to obtain the \(AC_s\), and was carried over when combined with the simulations used to obtain \(\hat{N}_{1,c}^M\).

3 Results and discussion

3.1 IPM assessment

The formulation of the capture-recapture part of the female and male IPMs was supported by the frequentist goodness-of-fit tests, which did not indicate lack of fit of Arnason-Schwarz multi-state model structures. Included in the IPM likelihoods, convergence in the MCMC chains was successfully achieved, with \(\hat{R}\) test values under 1.1. All the posterior estimates of demographic rates had a maximum proportion of the posterior distribution with the same sign as the mean (≥ 0.95). When fitted as independent capture-recapture likelihoods, the posterior distribution of survival and transition probability estimates showed maximum overlap compared to their distributions as part of the IPMs, which suggest very little discrepancy in their posterior distributions.

Model selection of female IPMs improved with correlated temporal random effects in \(\phi_B\), \(\phi_N\) and \(\psi_{BN}\), independent temporal random effects in survival and breeding probabilities (recruitment) of immature seals, and recapture probabilities modelled as invariant over time. Models without temporal effects on survival and breeding had much lower support. Due to limited sample size, the best male IPM had a simplified parameter structure, with temporal random effects for \(\phi_T\) and time-invariant options for the rest of the parameters. This structure was identical to the recent analyses of a reduce version of this data set and additional detail can be found in Forcada et al. (2023b).

The fit of GAMM likelihoods within the IPMs showed comparable support for models with temporal random effects and temporal autocorrelation (AR1), with marginal differences in WAIC. Estimated autocorrelation parameters were non-significantly different from 0, and thus a simplified Poisson-log-normal was selected throughout. When fitted independently, the GAMM likelihoods provided almost identical posterior distributions for parameters of interest as those obtained when fitted in the IPMs. This suggests that with a good choice of prior distribution and initial parameter estimates, fitting them independently was acceptable, with the advantage of saving computing time and fitting problems due to model complexity.

3.2 Life tables and generation length

The population projection matrix for females was slightly modified to include recruitment ages from three to seven, allowing for an absorbing state of seals aged seven or older never observed recruiting; this was a slightly different formulation than in the IPM, to derive life table and generation length estimates. The matrix models comprised averaged vital rates from 2001 to 2008, which is a pre-perturbation period before the major 2009 decline, as required by IUCN guidelines (IUCN 2024). The mean estimate and standard deviation are

\(\hat{A}_M=\begin{bmatrix} 0 &0 &0.021 &0.076 &0.129 &0.096 &0.085 &0.177 &0.263\\0.424 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0.974 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0.884 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0.435 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0.461 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0.442 &0.491 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0.309 &0.161\\0 &0 &0.053 &0.194 &0.327 &0.244 &0.216 &0.449 &0.667\end{bmatrix}\),

\(sd(\hat{A}_M)=\begin{bmatrix} 0 &0 &0.009 &0.012 &0.020 &0.026 &0.037 &0.018 &0.013\\0.035 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0.032 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0.055 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0.054 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0.066 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0.090 &0.140 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0.041 &0.011\\0 &0 &0.021 &0.029 &0.049 &0.065 &0.092 &0.040 &0.012\end{bmatrix}\).

The age structured models included a maximum observed age of 25, and the corresponding mean and standard deviation of the female life table is:

\(x\) \(px\) \(lx\) \(mx\) \(lxmx\) \(xlxmx\)
0 0.4237 (0.0357) 1 (0) 0 (0) 0 (0) 0 (0)
1 0.9744 (0.0312) 0.4237 (0.0357) 0 (0) 0 (0) 0 (0)
2 0.9369 (0.0527) 0.4125 (0.0329) 0.0206 (0.0091) 0.0085 (0.0037) 0.0169 (0.0075)
3 0.6396 (0.0601) 0.3860 (0.0328) 0.0871 (0.0127) 0.0335 (0.0050) 0.1005 (0.0149)
4 0.8007 (0.0480) 0.2462 (0.0257) 0.1759 (0.0170) 0.0432 (0.0054) 0.1728 (0.0214)
5 0.7693 (0.0400) 0.1969 (0.0213) 0.1950 (0.0163) 0.0383 (0.0046) 0.1916 (0.0230)
6 0.7917 (0.0328) 0.1512 (0.0162) 0.2130 (0.0158) 0.0322 (0.0039) 0.1931 (0.0235)
7 0.8023 (0.0183) 0.1197 (0.0131) 0.2234 (0.0157) 0.0267 (0.0033) 0.1870 (0.0234)
8 0.8073 (0.0129) 0.0961 (0.0113) 0.2295 (0.0157) 0.0220 (0.0028) 0.1762 (0.0228)
9 0.8097 (0.0107) 0.0776 (0.0097) 0.2332 (0.0156) 0.0181 (0.0024) 0.1627 (0.0218)
10 0.8108 (0.0097) 0.0629 (0.0083) 0.2355 (0.0153) 0.0148 (0.0021) 0.1479 (0.0206)
11 0.8113 (0.0092) 0.0510 (0.0071) 0.2370 (0.0151) 0.0121 (0.0017) 0.1329 (0.0192)
12 0.8116 (0.0089) 0.0414 (0.0061) 0.2380 (0.0149) 0.0099 (0.0015) 0.1183 (0.0179)
13 0.8117 (0.0088) 0.0337 (0.0052) 0.2388 (0.0148) 0.0080 (0.0013) 0.1044 (0.0165)
14 0.8117 (0.0087) 0.0274 (0.0044) 0.2393 (0.0146) 0.0065 (0.0011) 0.0916 (0.0151)
15 0.8117 (0.0086) 0.0222 (0.0037) 0.2397 (0.0145) 0.0053 (0.0009) 0.0798 (0.0138)
16 0.8117 (0.0086) 0.0181 (0.0032) 0.2400 (0.0144) 0.0043 (0.0008) 0.0693 (0.0125)
17 0.8116 (0.0085) 0.0147 (0.0027) 0.2402 (0.0144) 0.0035 (0.0007) 0.0599 (0.0113)
18 0.8116 (0.0085) 0.0119 (0.0023) 0.2404 (0.0143) 0.0029 (0.0006) 0.0516 (0.0101)
19 0.8116 (0.0085) 0.0097 (0.0020) 0.2406 (0.0142) 0.0023 (0.0005) 0.0443 (0.0091)
20 0.8115 (0.0085) 0.0079 (0.0017) 0.2407 (0.0142) 0.0019 (0.0004) 0.0379 (0.0081)
21 0.8115 (0.0085) 0.0064 (0.0014) 0.2408 (0.0142) 0.0015 (0.0003) 0.0324 (0.0072)
22 0.8115 (0.0085) 0.0052 (0.0012) 0.2409 (0.0141) 0.0013 (0.0003) 0.0276 (0.0064)
23 0.8114 (0.0085) 0.0042 (0.0010) 0.2410 (0.0141) 0.0010 (0.0002) 0.0235 (0.0056)
24 0.8114 (0.0084) 0.0034 (0.0009) 0.2410 (0.0141) 0.0008 (0.0002) 0.0199 (0.0050)

Table 2. Antarctic fur seal female life table for pre-disturbance (pre-decline) years based on derived vital rate estimates from the fitted female IPM model. Column parameters are referred to interval \([x, x+1]\), where \(x\) is the age at the start of the age interval. These are \(px\), or probability of survival for seals entering the interval; \(lx\), or survivorship, which is the proportion of initial cohort surviving to the start of age interval; \(mx\), or per-capita rate of reproduction during interval; \(lxmx\), or expected number of offspring per cohort member produced during interval; and \(xlxmx\). Standard errors are in parentheses.


The population projection matrix for males was slightly modified to include recruitment ages from seven to eleven, allowing for an absorbing state of animals eleven or older never observed recruiting. It comprised averaged vital rates from 1995 to 2007, which is the pre-perturbation period before the major 2009 decline. For this period the male population growth was near stable (Forcada et al. 2023b). The mean estimate and standard deviation of the matrix model are


\(\hat{A}_M=\begin{bmatrix} 0 &0 &0 &0 &0 &0 &0.003 &0.017 &0.038 &0.052 &0.094 &0.136 &0.182\\0.719 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0.813 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0.840 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0.83 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0.790 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0.721 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0.620 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0.488 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0 &0.358 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0 &0 &0.295 &0.212 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0.171 &0.142\\0 &0 &0 &0 &0 &0 &0.009 &0.047 &0.108 &0.149 &0.267 &0.387 &0.518\end{bmatrix}\),

\(sd(\hat{A}_M)=\begin{bmatrix} 0 &0 &0 &0 &0 &0 &0.002 &0.005 &0.012 &0.023 &0.053 &0.022 &0.021\\0.014 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0.013 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0.010 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0.010 &0 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0.016 &0 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0.024 &0 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0.032 &0 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0.049 &0 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0 &0.079 &0 &0 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0 &0 &0.115 &0.123 &0 &0\\0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0 &0.057 &0.053\\0 &0 &0 &0 &0 &0 &0.006 &0.014 &0.032 &0.066 &0.149 &0.061 &0.056\end{bmatrix}\).

This model includes a reproductive component (first row), which is required to obtain a life table for estimating generation length based on the “cohort” method of Coale (1972). The mean and standard deviation of the derived life table, including a maximum observed age of 18, is

\(x\) \(px\) \(lx\) \(mx\) \(lxmx\) \(xlxmx\)
0 0.7193 (0.0145) 1 (0) 0 (0) 0 (0) 0 (0)
1 0.8129 (0.0131) 0.7193 (0.0145) 0 (0) 0 (0) 0 (0)
2 0.8401 (0.0099) 0.5849 (0.0210) 0 (0) 0 (0) 0 (0)
3 0.8299 (0.0106) 0.4915 (0.0226) 0 (0) 0 (0) 0 (0)
4 0.7890 (0.0167) 0.4080 (0.0210) 0 (0) 0 (0) 0 (0)
5 0.7190 (0.0243) 0.3219 (0.0180) 0 (0) 0 (0) 0 (0)
6 0.6271 (0.0312) 0.2315 (0.0158) 0.0033 (0.0021) 0.0008 (0.0005) 0.0045 (0.0029)
7 0.5353 (0.0501) 0.1453 (0.0141) 0.0192 (0.0048) 0.0028 (0.0008) 0.0196 (0.0056)
8 0.4859 (0.0870) 0.0781 (0.0121) 0.0528 (0.0107) 0.0042 (0.0012) 0.0333 (0.0099)
9 0.5201 (0.1058) 0.0385 (0.0113) 0.0962 (0.0160) 0.0038 (0.0015) 0.0341 (0.0138)
10 0.6042 (0.0813) 0.0210 (0.0099) 0.1498 (0.0215) 0.0032 (0.0018) 0.0323 (0.0183)
11 0.6377 (0.0358) 0.0133 (0.0081) 0.1665 (0.0199) 0.0022 (0.0014) 0.0246 (0.0158)
12 0.6409 (0.0252) 0.0087 (0.0059) 0.1702 (0.0202) 0.0015 (0.0010) 0.0179 (0.0124)
13 0.6384 (0.0414) 0.0057 (0.0041) 0.1713 (0.0204) 0.0010 (0.0007) 0.0127 (0.0093)
14 0.6143 (0.1273) 0.0037 (0.0028) 0.1718 (0.0205) 0.0006 (0.0005) 0.0089 (0.0068)
15 0.5454 (0.2295) 0.0025 (0.0019) 0.1721 (0.0207) 0.0004 (0.0003) 0.0064 (0.0049)
16 0.4597 (0.2909) 0.0018 (0.0013) 0.1731 (0.0207) 0.0003 (0.0002) 0.0049 (0.0036)
17 0.3759 (0.3189) 0.0014 (0.0009) 0.1742 (0.0212) 0.0002 (0.0002) 0.0042 (0.0026)
18 0 (0) 0.0012 (0.0006) 0.1759 (0.0214) 0.0002 (0.0001) 0.0037 (0.0020)

Table 3. Antarctic fur seal male life table for pre-disturbance (pre-decline) years based on derived vital rate estimates from the fitted female IPM model. Column parameters are referred to interval \([x, x+1]\), where \(x\) is the age at the start of the age interval. These are \(px\), or probability of survival for seals entering the interval; \(lx\), or survivorship, which is the proportion of initial cohort surviving to the start of age interval; \(mx\), or per-capita rate of reproduction during interval; \(lxmx\), or expected number of offspring per cohort member produced during interval; and \(xlxmx\). Standard errors are in parentheses.

From the sex-specific life tables, the cohort-method generation length estimate is 8.21 (SD=0.28) for females and 9.90 (SD=0.58) for males. Using mean absolute abundance estimates for mature females and mature males between 2007-2009, when the population of South Georgia was estimated to be at its maximum, the mean generation length weighted by the abundance of mature seals of each sex is 8.63 (SD=0.26). This value is used subsequently for the population reduction assessment.

3.3 Adjustments to scale-up pup production to mature and total population

Using the IPM estimates of absolute abundance, an estimate of the ratio of mature individuals to pup production for the interval 2001 and 2005 was 1.482 [95% CI: 1.392, 1.582], and was used as the mean adjustment from a pup count to mature population. The ratio of total population, including mature, pre-breeders and pre-territorials, to pup production, equivalent to a mean adjustment to total population, was 3.026 [2.767, 3.303].

Previously, Payne (1979) indicated that the total population size at the end of the breeding season was 4.1 times the number of pups born. However, to account for unobservable population components, including, for example, the majority of the pre-territorial males which is nearly 80% of the total male population (Forcada et al. 2023b), he could only use limited data. His adjustment was obtained from calculations of body mass at age from a small sample of known age seals from Bird Island, with the remark: “Nothing is known of the structure of the male herd other than the mean age of a small sample of breeding bulls” (Payne 1979). He also used for his adjustment a population growth rate calculated from pup production values obtained over a decade apart, between late 1950s and mid 1970s, with substantially different and hardly comparable methods (Croxall and Prince 1979). Additionally, he assumed survival rates at age, for a population age structure obtained with little data that was largely affected by unquantified loss of plastic tags, which were the main individual marking method at the time. As a result, his highly biased multiplier is not even within the confidence limits of the recently and robustly estimated factor of 2.96 [2.64, 3.33] for South Georgia (Forcada et al. 2023b); although it has been repeatedly used until recently (e.g., Boyd 1993; Hofmeyr et al. 2005; Krause et al. 2022). For this analysis, I updated pup count adjustments with the IPM models that use the most up-to-date recapture data of known females (n=4,029) and males (n=6,108), and all of them were reliably tagged with electronic chips and/or genetic markers, which are permanent marking methods.

3.6 Global and subpopulation assessment using Red List category and criteria

3.6.1 South Georgia and South Orkney Islands

Using the IPM simulations, I derived an estimated mean annual change from 1996 to 2006, which was 0.98 (SD=0.01). It was consistent with a previously estimated stochastic \(log\lambda\) for the same period at Bird Island (Forcada, Trathan, and Murphy 2008), and therefore, I used it to project the population of South Georgia from 1999 (\(t_{3G}\)) to 2007/9 (\(t_1\)). For the observed period, \(t_1\) to \(t_T\) (2022), the annual rate of decline was faster, -3.93% [CI: -2.73; -5.13], and the estimated reduction over three generations was -0.57 [-0.65; -0.48], with a probability of 0.92 of being lower than or equal to -0.50. This supported a higher than 50% reduction, and the Red List Category Endangered as the most likely, based on the A2 criterion assessment, as discussed above.

At the South Orkney Islands, pup counts over the last three decades have been zero or near zero at the long-term studies of Laurie Island and Signy Island (Carlini et al. 2006; Waluda, Gregory, and Dunn 2010); and have been zero at Signy Island since 2010. Surveys of Gosling, Monroe, Powell, and Michelsen Islands in 1974 produced a total count of 65 pups (Laws 1981), but data for these areas are lacking because they are seldom visited for seal surveys. Based on the current absence of puping at study sites, my assumption is that fur seal pupping at the South Orkney Island group continues to be marginal, as discussed in Dunn et al. (2025).

3.6.2 South Sandwich Islands

Data from surveys of the South Sandwich Islands archipelago conducted in January and February 1997, and completed during February of 1998 (Convey, Morton, and Poncet 1999), provide the best baseline for total pup counts at this difficult to access location. The counts did not account for mortality and were grossly rounded; but better, or any subsequent data for that matter, are unavailable for most of the islands. Most of the pups in 1997/8 were found in Zavodovski (n=500), Visokoi (n=900), and Candlemas Islands (n=346), at the northwest of the Archipelago. In December 2023, a fur seal survey of Zavodovski Island only was completed, with an estimated RPAS count of 1,265 pups.

To assess the staus of this breeding location, and given the observed increase at Zavodovski Island, I assumed no change at other unsurveyed locations. Thus, the mature population in \(t_{3G}\) from a total count of 1,750 pups in 1997/8 was 2,535 (SD=93); and in 2025 was estimated as 3,782 (SD=124). This implies a category of Least concern, on the basis of the population increase at Zavodovski Island only, but with two major caveats. First, 1997/1998 saw the second strongest El Niño event on record, which brought about significant breeding deferral and low pup production at South Georgia (Forcada et al. 2005; Forcada, Trathan, and Murphy 2008). This was likely due to direct climatic effects on the marine environment around nearby South Georgia and adjacent areas (Murphy et al. 2007). In all certainty, it would also be associated with a lower than average pup count at the South Sandwich Islands (Forcada et al. 2005). And second, the assessment is questionable due to lack of data of at least Visokoi Island, where pups were most abundant in 1998. So for the South Sandwich Islands, the most realistic category is Data Deficient.

3.6.3 South Shetland Islands

Using a total pup count of 9,530 (SD=135) of the South Shetland Islands archipelago obtained in 1995/6 (Meyer, Walker, and Holt 1996), and a total count of 10,057 (SD=142) from 2001/2 (\(t_1\)) (Table 1), I estimated an annual growth rate of 1.009 (SD=0.003) for the period 1996-2002. I then used this rate to back-project the population from \(t_1\) to \(t_{3G}\), and obtained a mature population size in \(t_{3G}\) (1999) of 14,531 (SD=556). In \(t_1\) (2002), the mature popuation size was estimated as 14,913 (SD=531), and as 1,336 (SD=44) in \(t_T\) (2025). The annual rate of change for the observed period (2002 to 2025) was -9.96% [-9.85; -10.06], and the estimated reduction over three generations was -0.908 [-0.911; -0.904], with maximum probability of being higher than 80%, and equal or higher than 90%. This corresponds to a category of Critically Endangered under Criterion A2, and also under Criterion A1, as the reduction was higher than 90%; although the causes of the decline are only partially known, and have not ceased and are unlikely to be reversible, so Criterion A2 applies.

3.6.4 Bouvetøya

The pup count at Nyrøysa (Bouvetøya) for 1999, which is both \(t_1\) and \(t_{3G}\), was 15,448 (SD=370) (Hofmeyr et al. 2005), and corresponded to a mature population size of 22,832 (SD=1,141). A pup production estimate of 8,080 (SD=310) during season 2015 (\(t_T\)) was the latest data available, as research teams have not returned to the island since (A. Lowther, personal communication). This corresponds to a mature population size of 11,990 (SD=610), and an annual rate of change of -3.98% [CI:3.44; 4.52] between 1999 and 2015. Assuming a similar rate of decline until the assessment year (\(t_{AY}\)=2025), the reduction over three generations is -0.65 [-0.70; -0.60], which corresponds with maximum probability to a category of Endangered under criterion A2.

3.6.5 East Antarctica

3.6.5.1 Marion Island and Prince Edward Island

For Marion Island, the closest available pup count to 1999 (\(t_{3G}\)) was obtained in 2000 (\(t_1\)) and was 343 (Wege et al. 2016), which is equivalent to an estimated mature population of 509 (SD=17), and a mature population in \(t_{3G}\) of 470 (SD=36). The latest pup production estimate is for 2025 (\(t_T\)), of 3,172 (SD=156) pups (Table 1), and corresponds to a mature population of 4,704 (SD=278). The estimated annual increase is 9.30% [8.88; 9.72], which supports a Red List Category of Least Concern.

There are no data for Prince Edward Island since 2001, when a pup count of 404 was reported (Bester, Ryan, and M. 2003). This corresponded to an estimated mature population of 599 [CL: 563; 639]. Without any data to evaluate \(N^M_{AY}\), I used a conservative approach to contribute to the evaluation of the mature population of the East Antarctica subpopulation, and assumed no population change at this breeding location. I used this approach for other data deficient locations (see below).

3.6.5.2 Crozet Islands

The main breeding location at Crozet Islands archipelago is in Possession Island, and has been surveyed repeatedly (Guinet et al. 2025). The pup count in 2000 was 234, which corresponds to a mature population of 347 (SD=11), and a back-projected mature population of 337 (SD=15) in \(t_{3G}\). The latest pup count was in 2024, of 530 pups, and corresponded to a mature population of 786 (SD=26). The projected population in \(t_{AY}\) is 813 (SD=27), assuming the same estimated annual growth of 3.47% obtained for the observed period. This supports a Red List Category of Least Concern.

3.6.5.3 Kerguelen Islands

There are several breeding locations in the Kerguelen Islands archipelago, but only the Courbet Peninsula has been regularly counted (Guinet et al. 2025). The closest assessment to 1999 (\(t_{3G}\)) was in 2000, of 1,600 pups. There is another rounded pup count of 3,600 pups in 2000 at Nuageuses Islands (Guinet et al. 2025), where the numbers are likely to have remained stable (C. Guinet, personal communication); but there are no more recent counts. Assuming that the population at Nuageuses Islands has remained constant, I added the pup count of 3,600 to both of two pup counts obtained at Courbet Peninsula to facilitate the assessment of the entire East Antarctica subpopulation. The mature population in \(t_{3G}\) would then be equivalent to 7,616 (SD=267), and with the latest pup count at Possession Islands of 4,793, it would be 12,698 (SD=415) in \(t_{AY}\). The resulting estimated annual increase is 2.02%, which supports a Red List Category of Least Concern.

3.6.5.4 Heard Island and McDonald Islands

Due to its remoteness and difficult access, the Heard Island and McDonald Islands group have been seldom visited, let alone counted the fur seal pups being born on them. In March 1980, “up to a 100 pups and 200 other” (ditto, Johnstone (1982)) were counted at McDonald Islands, being the last available data for this assessment (Johnstone 1982); although the population was thought to be increasing (SCAR 2006). As a result, this breeding location continues to be Data Deficient, and because the data preceds \(t_{3G}\) by far, it is not considered numerically in the assessment of East Antarctica.

In contrast, Heard Island was visited on several occasions, the more relevant being 2001, when 1,012 pups were counted, and the last one in 2004, with a count of 1,278 pups (S. Goldsworthy, personal communication). This would correspond to a mature population of 1,501 (SD=49) in 2001, and 1,895 (SD=62) in 2004, with a possible population in \(t_3G\) of 1,296 (SD=90). For the observed period, the annual increase was 8.089%. However, these data alone cannot warrant a new assessment, and I used the same approach as MacDonald and Kerguelen Islands, keeping the last value (in this case for 2004) stable until \(t_{AY}\). The Red List category for Heard Island is therefore Data Deficient.

3.6.5.5 Macquarie Island

There are three species of fur seals that breed at Macquarie Island; Antarctic, subantarctic (A. tropicalis) and New Zealand (A. forsteri) fur seals, although there are no females of New Zealand fur seals. All species hybridise, and best attempts at discriminating pups to species level depend on genetic rather than phenotypic identification, which tends to overestimate A. gazella as species (Goldsworthy et al. 2009). Based on these early counts, the number of A. gazella pups were approximately 100-105 in 1999, which would give an \(N^M_{3G}\) of 156 (SD=498), and an \(N^M_T\) of 280 (SD=9) in 2012, when long-term studies were terminated. Annual change between these estimates and observed period would be 4.65%. As with previous populations, this supports a Data Deficient category; and I used the assumption of stability after the last count to evaluate its contribution to the East Antarctica subpopulation.

Breeding location \(N^M_{3G}\) \(N^M_{t_1}\) \(N^M_{t_T}\) \(N^M_{AY}\) Annual change Category
Marion Island 470 (36) 4,704 (278) 9.30% (0.002) Least Concern
Prince Edward Is. 2001; 599 (20) Data Deficient
Crozet Islands 337 (15) 2001: 347 (11) 2024: 786 (26) 813 (27) 3.46% Least Concern
Kerguelen Islands 7,616 (267) 2000: 7,712 (252) 2024: 12,448 (407) 12,698 (415) 2.02% Least Concern
Heard Island 1,296 (90) 2001: 1,501 (49) 2004: 1,895 (62) Data Deficient
Macquarie Island 156 (8) 2012: 280 (9) Data Deficient

Table 6. Population assessment and IUCN Red List Category at breeding locations of the East Antarctica subpopulation. Column \(N^M_{3G}\) is population size in 1999, three times the generation length or 26 years; \(N^M_{t_1}\) is the first available population size; \(N^M_{t_T}\) is the second available population size; and \(N^M_{AY}\) is the current population size. Values in parentheses are standard deviations derived from Bayesian simulations. Cells in columns for \(N^M_{t_1}\) and \(N^M_{t_T}\) start with the year when the value was obtained.

3.6.5.6 East Antarctica subpopulation

For the assessment of East Antarctica (Table 6), I used the sum of the estimates of \(N^M_{3G}\) for all breeding locations of this subpopulation except McDonald Islands. To obtain \(N^M_{AY}\), I used the estimates from the most recent data, or their projection to \(AY\) if these two values were not the same; and for other cases, I maintained the value of the more recent of \(N^M_{t_1}\) or \(N^M_{t_T}\) to provide an \(N^M_{AY}\), assuming stability in growth or no change until the assessment year. This was the case of Maquarie Island, Heard Island, part of Kerguelen Islands, and Prince Edward Islands.

As seen above, when data were available for \(t_{3G}\) and \(t_{AY}\), population growth was positive; and all the locations which had two population estimates, including those which were Data Deficient, also showed positive growth. Locations that were Data Deficient contributed little to total pup production/mature population, except perhaps Heard Island. This suggests that keeping their estimates constant until \(t_{AY}\) was likely to have a minor effect in the assessment of the entire subpopulation, with the Red List category unchanged; i.e., Least Concern. The annual combined increase for East Antarctica was 2.62% [CI: 2.45; 2.80] (Table 7).

3.6.6 Global population and species

From the combined estimates of \(N^M_{3G}\) and \(N^M_{AY}\) of all subpopulations, including South Sandwich Islands, the annual change for the global population was estimated as -3.19% [CI:-4.015; -2.347]. The estimated reduction was -0.57 [-0.64; -0.48], with a probability of 0.9 of being equal or lower than -0.50 (50%); which corresponds to to a Red List Category of Endangered (Table 7). This implies a major change of the conservation status of the species, which was previously considered of Least Concern but now is in an extinction risk Category.

Currently, only the East Antarctica subpopulation continues to be Least Concern, with the caveat that some of the breeding locations, although not generally the most productive, can be considered Data Deficient. In contrast, the South Atlantic region has seen the most precipitous declines for the species in the last decades. Numerically, South Georgia has had an overwhelming influence both in the reduction of the global population and also the population of this region. This is because it comprises approximately 96% of the global population according to this new assessment (Table 7).

In terms of severity of the reduction, the South Shetland Islands is the most affected subpopulation, as it can now be considered Critically Endangered, and its size is now smaller than at some of the East Antarctica breeding locations; whereas before it was the third largest subpopulation. It is potentially smaller in size than the population at the South Sandwich Islands, but concerns on data availability there preclude a proper comparison. The genetic identity of the South Sandwich Islands remains to be examined in order to understand if it can be considered a different subpopulation; and a new complete population assessment is needed.

To put this into perspective, the South Shetland Islands subpopulation which has been considered a genetic reservoir (Krause et al. 2024), is more likely to be now at risk of disappearing, or even of introgression through mixing with other subpopulations, given a potential temporary migratory expansion southwards of males and females (see Larsen et al. 2025) as the habitat changes (Ouled-Cheikh et al. 2024). This will provide different options for dispersal, recolonisation, and potentially new breeding locations, although this may only occur if the appearance of new habitats further South with sea ice loss also brings about appropriate conditions to accommodate the stringent seasonal central place foraging demands of breeding females.

Population \(N^M_{3G}\) \(N^M_{AY}\) % of global Population change Red List Category
Annual change Reduction (R) P(R ≤ -0.5)
South Georgia 2,136,953 (153,084) 910,114 (89,753) 96.40 (0.37) -3.24% (0.44) -0.57 (0.05) R ≤ -0.5: 0.92 Endangered
South Sandwich Islands 474 (34) 3,782 (124)* 0.40 (0.04) 1.40% Data Deficient
South Shetland Islands 14,529 (551) 1,336 (44) 0.14 (0.01) -8.77% (0.07) -0.908 (0.002) R ≤ -0.5: 1 R ≤ -0.8: 1 Critically Endangered
Bouvetøya 22,832 (1,141) 8,009 (582) 0.86 (0.11) -3.95% (0.29) -0.65 (0.03) R ≤ -0.5: 1 Endangered
East Antarctica 10,474 (410) 20,514 (710) 2.19 (0.22) 2.62% (0.09) Least Concern
GLOBAL 2,187,424 (153,141) 943,756 (89,727) 100 -3.19% (0.43) -0.57 (0.05) R ≤ -0.5: 0.90 Endangered

Table 7. Population assessment and IUCN Red List Category for Antarctic fur seal global population and subpopulations. Columns \(N^M_{3G}\) and \(N^M_{AY}\)) are population estimates for the assessment period (1999-2025), equivalent to three times the generation length (3G) which is 26 years for the assessment year of 2025. The % of global value is the percentage of the size of each subpopulation, or breeding location in the case of the South Sandwich Islands, with respect to the global population. The probability of a population size reduction higher than 50% (or 80%) is under P(R ≤-0.5). Values in parentheses are standard deviations derived from Bayesian simulations. *The \(N^M_{AY}\) value for South Sandwich Islands is based on insufficient data; and the genetic identity of the the South Sandwich Islands is unknown, and for this assessment is assumed to not be the same as the South Georgia population.

3.7 Justification of the change in status

The results of this new assessment support a status of Endangered according to criterion A2 for the Antarctic fur seal, which was previously listed as least concern. Based on current studies, the causes of the reduction are highly unlikely to be reversible (e.g., due to insufficient climate change mitigation), have not ceased, and some are not yet understood. As noted above, the ongoing declines largely respond to directional environmental effects, particularly at South Georgia (Forcada et al. 2023b), where increasing temperature and cascading ecosystem-wide repercussions determines the staple food supply required to sustain such a vast fur seal population (Atkinson et al. 2019; Hill, Phillips, and Atkinson 2013; Kawaguchi et al. 2024; McBride et al. 2021). But the detail of the mechansims and processes remains elusive. In addition, there are competing effects such as local top-down control at the South Shetland Islands (Krause et al. 2022), and potential effects for which quantification is almost intractable due to lack of data (see Introduction).

The large population reduction estimated for South Georgia has an overwhelming effect on the species Category, due to the numerical importance of this subpopulation. For long-lived species, reductions for criterion A are calculated over 3 generations, to avoid confusion with fluctuations, but data from shorter time series can be used to calculate the 3-generation reduction (IUCN 2024). Here, the time series data for the projection of South Georgia spans 2 generations (2007-2022), but there is very strong supporting evidence of the long-term decline at South Georgia. This is shown by data from Bird Island, which has been declining since 1990 (Forcada, Trathan, and Murphy 2008) (see also Fig. 4, central panel), a period of over 4 generations. The long-term decline at South Georgia is also supported by the population index trajectory of Signy Island (South Orkney Islands), which is driven by fur seals from South Georgia in its majority (Boyd et al. 1998), and its variation is explained by population changes at Bird Island (Fig. 5). Signy Island follows the same trajectory as Bird Island when mid-term sea ice effects are taken into account (Dunn et al. 2025). Here, the time-series illustrates well an example of fluctuation due to environmental effects (2000-2009), embedded in an ongoing long-term decline (1994-2025) (Fig.4, lower panel; Fig. 6).

This assessment supports an important change, due entirely to newly obtained data, which supersedes previously used insufficient and incorrect data. The newly proposed Red List status can be thus considered a Nongenuine change, according to the latest IUCN definition of change [page 13; IUCN (2024)]. “New information - The change in category is the result of better knowledge about the taxon, e.g. owing to new or newly synthesized information about the status of the taxon (e.g., better estimates for population size, range size or rate of decline)”, which in this case is better estimates for population size, and rate of decline. And “Incorrect data - The previous category was applied in error because incorrect data were used”. In this case, previously made incorrect assumptions on population size and growth, particularly for South Georgia.

The previous assessment was based on incomplete data which precluded an investigation of population reduction. Survey data, particularly from South Georgia which comprises 96% of the global abundance, were unavailable at the time of the previous assessment. Instead, the assessment quoted old and non-validated abundance estimates for South Georgia for 1991 (SCAR 2006), which were very different from those reported previously (Boyd 1993) (2.7 million compared to 1.5 million, respectively); and from which population values for 1999/00 were derived as 4.5 and 6.2 million seals (SCAR 2006). These figures were allegedly obtained from female counts, but there are no records or any other evidence of any female counts or surveys of South Georgia in 1999/00, at the British Antarctic Survey or elsewhere. Recently, the population in 2007/9, at the maximum of its recovery at South Georgia, has been robustly estimated as less than 4 million (3.5 million [CI:3,1 - 3,9]). It is likely that the previous incorrect estimates were extrapolated from female counts in 1991, reevaluated with some adjustment factor, and assuming an implausible annual growth rate of 6% to 14% (SCAR 2006) obtained at Bird Island, the only place were long-term data were available. As seen above, this growth rate is inconsistent with the much lower annual growth rates observed at Bird Island [see also this analysis; Forcada et al. (2023b)].

Therefore, without representative pup counts or other survey data, previously reported numbers were conjectural uneducated guesses, only supported by a personal communication in SCAR (2006). This, as dicussed in Forcada et al. (2023b), led to a category of Least Concern in 2016 (Hofmeyr 2016), where the mature female population was reported to be around 550,000, a value which was not in the references cited to obtain such estimate. Additionally, this population value was cited as questionable by a personal opinion from a blog which was, again, based on no data and pure speculative criticism, adding even more confusion. For this reason, the previous assessment is not based on reliable population data, and therefore, the new assessment reports a Nongenuine change, due to the use of new and robust assessments of the latest data.

4 References

Abadi, F., O. Gimenez, R. Arlettaz, and M. Schaub. 2010. “An Assessment of Integrated Population Models: Bias, Accuracy, and Violation of the Assumption of Independence.” Ecology 91 (1): 7–14. https://doi.org/https://doi.org/10.1890/08-2235.1.
Agnew, D. J. 1997. “The CCAMLR Ecosystem Monitoring Programme.” Journal Article. Antarctic Science 9 (3): 235–42. <Go to ISI>://WOS:A1997XV76100002.
Atkinson, Angus, Simeon L. Hill, Evgeny A. Pakhomov, Volker Siegel, Christian S. Reiss, Valerie J. Loeb, Deborah K. Steinberg, et al. 2019. “Krill (Euphausia Superba) Distribution Contracts Southward During Rapid Regional Warming.” Journal Article. Nature Climate Change. https://doi.org/10.1038/s41558-018-0370-z.
Aznar-Alemany, Òscar, Xuefei Yang, Mariana B. Alonso, Erli Schneider Costa, João Paulo M. Torres, Olaf Malm, Damià Barceló, and Ethel Eljarrat. 2019. “Preliminary Study of Long-Range Transport of Halogenated Flame Retardants Using Antarctic Marine Mammals.” Science of The Total Environment 650: 1889–97. https://doi.org/https://doi.org/10.1016/j.scitotenv.2018.09.287.
Baines, Mick, Jennifer A. Jackson, Sophie Fielding, Vicky Warwick-Evans, Maren Reichelt, Claire Lacey, Simon Pinder, and Philip N. Trathan. 2022. “Ecological Interactions Between Antarctic Krill (Euphausia Superba) and Baleen Whales in the South Sandwich Islands Region – Exploring Predator-Prey Biomass Ratios.” Deep Sea Research Part I: Oceanographic Research Papers 189: 103867. https://doi.org/https://doi.org/10.1016/j.dsr.2022.103867.
Bennison, A., A. Stacey, Banyard A, F. Blockley, M. Blyth, E. Browne, G. Day, et al. 2024. “A Case Study of Highly Pathogenic Avian Influenza (HPAI) H5N1 at Bird Island, South Georgia: The First Documented Outbreak in the Subantarctic Region.” Journal Article. Bird Study 71 (4): 380–91. https://doi.org/10.1080/00063657.2024.2396563.
Besbeas, P., S. N. Freeman, B. J. T. Morgan, and E. A. Catchpole. 2002. “Integrating Mark–Recapture–Recovery and Census Data to Estimate Animal Abundance and Demographic Parameters.” Biometrics 58 (3): 540–47. https://doi.org/https://doi.org/10.1111/j.0006-341X.2002.00540.x.
Bester, M. N., P. G. Ryan, and Dyer B. M. 2003. “Population Numbers of Fur Seals at Prince Edward Island, Southern Ocean.” African Journal of Marine Science 25 (1): 549–54. https://doi.org/10.2989/18142320309504044.
Biuw, Martin, Ulf Lindstrøm, Jennifer A. Jackson, Mick Baines, Nat Kelly, George McCallum, Georg Skaret, and Bjørn A. Krafft. 2024. “Estimated Summer Abundance and Krill Consumption of Fin Whales Throughout the Scotia Sea During the 2018/2019 Summer Season.” Journal Article. Scientific Reports 14 (1): 7493. https://doi.org/10.1038/s41598-024-57378-3.
Bonner, W. N. 1968. “The Fur Seal of South Georgia.” Report.
Boyd, I. L. 1993. “Pup Production and Distribution of Breeding Antarctic Fur Seals (Arctocephalus Gazella) at South Georgia.” Antarctic Science 5 (1): 17–24. https://doi.org/10.1017/S0954102093000045.
Boyd, I. L., D. J. McCafferty, K. Reid, R. Taylor, and T. R. Walker. 1998. “Dispersal of Male and Female Antarctic Fur Seals (Arctocephalus Gazella).” Canadian Journal Fisheries Aquatic Sciences 55 (4): 845–52.
Brault, Emily K., Michael E. Goebel, Heidi N. Geisz, Elizabeth A. Canuel, and Rebecca M. Dickhut. 2013. “Inter-Annual Variation of Persistent Organic Pollutants (POPS) in an Antarctic Top Predator Arctocephalus Gazella.” Journal Article. Environmental Science & Technology 47 (22): 12744–52. https://doi.org/10.1021/es403577f.
Brownell, R. L., D. J. Krause, A. M. M. Baylis, C. A. Bonin, L. R. Oliveira, M. M. Uhart, M. Ulloa, and G .M. Watters. 2024. “Avian Influenza H5N1 Threatens Imperiled Krill-Dependent Predators in Antarctica.” Journal Article. Frontiers in Marine Science Volume 11 - 2024. https://doi.org/10.3389/fmars.2024.1453737.
Carlini, A. R., G. A. Daneri, R. Casaux, and M. E. I. Márquez. 2006. “Haul-Out Pattern of Itinerant Male Antarctic Fur Seals (Arctocephalus Gazella) at Laurie Island, South Orkney Islands.” Polar Research 25 (2): 139–44. https://doi.org/10.3402/polar.v25i2.6243.
Caswell, H. 2001. Matrix Population Models : Construction, Analysis, and Interpretation. Sinauer Associates. Book. 2nd ed. Sunderland, Mass.: Sinauer Associates.
Cazelles, B., M. Chavez, D. Berteaux, F. Menard, J. O. Vik, S. Jenouvrier, and N. C. Stenseth. 2008. “Wavelet Analysis of Ecological Time Series.” Journal Article. Oecologia 156 (2): 287–304. https://doi.org/10.1007/s00442-008-0993-2.
Coale, A. J. 1972. The Growth and Structure of Human Populations: A Mathematical Investigation. Book. Princeton University Press. https://books.google.co.uk/books?id=PlaZQgAACAAJ.
Convey, P., A. Morton, and J. Poncet. 1999. “Survey of Marine Birds and Mammals of the South Sandwich Islands.” Polar Record 35 (193): 107–24. https://doi.org/10.1017/S0032247400026450.
Costa, D. P., M. J. Weise, and J. P. Y. Arnould. 2007. “Potential Influences of Whaling on the Status and Trends of Pinniped Populations.” In Whales, Whaling, and Ocean Ecosystems. University of California Press. https://doi.org/10.1525/california/9780520248847.003.0027.
Croxall, J. P., T. S. McCann, P. A. Prince, and P. Rothery. 1988. “Reproductive Performance of Seabirds and Seals at South Georgia and Signy Island, South Orkney Islands, 1976–1987: Implications for Southern Ocean Monitoring Studies.” Conference Proceedings. In, 261–85. Antarctic Ocean and Resources Variability. Springer Berlin Heidelberg.
Croxall, J. P., and P. A. Prince. 1979. “Antarctic Seabird and Seal Monitoring Studies.” Polar Record 19 (123): 573–95. https://doi.org/10.1017/S0032247400002680.
Davis, R. A., and R. Wu. 2009. “A Negative Binomial Model for Time Series of Counts.” Biometrika 96 (3): 735–49. https://doi.org/10.1093/biomet/asp029.
Dennis, B., and M. L. Taper. 1994. “Density Dependence in Time Series Observations of Natural Populations: Estimation and Testing.” Ecological Monographs 64 (2): 205–24. https://doi.org/https://doi.org/10.2307/2937041.
Doidge, D. W., J. P. Croxall, and J. R. Baker. 1984. “Density-Dependent Pup Mortality in the Antarctic Fur Seal Arctocephalus Gazella at South Georgia.” Journal of Zoology 202 (3): 449–60. https://doi.org/https://doi.org/10.1111/j.1469-7998.1984.tb05095.x.
Duck, C. D. 1990. “Annual Variation in the Timing of Reproduction in Antarctic Fur Seals, Arctocephalus Gazella, at Bird Island, South Georgia.” Journal of Zoology 222 (1): 103–16. https://doi.org/https://doi.org/10.1111/j.1469-7998.1990.tb04032.x.
Dunn, M. J., C. M. Waluda, S. Adlard, D. Fox, A. S. Lynnes, T. I. Morley, and J. Forcada. 2025. “Temporary Absence of Warming in the Northern Weddell Sea Validates Expected Responses of Antarctic Seals to Sea Ice Change.” Global Change Biology 31 (6): e70290. https://doi.org/https://doi.org/10.1111/gcb.70290.
Forcada, J. 2021. “Antarctic Fur Seal (Arctocephalus Gazella, Peters 1875).” In Ecology and Conservation of Pinnipeds in Latin America, edited by G. Heckel and Y. Schramm, 51–76. Cham: Springer International Publishing. https://doi.org/10.1007/978-3-030-63177-2_4.
Forcada, J., and J. I. Hoffman. 2014. “Climate Change Selects for Heterozygosity in a Declining Fur Seal Population.” Nature 511 (7510): 462–65. https://doi.org/10.1038/nature13542.
Forcada, J., J. I. Hoffman, O. Gimenez, I. J. Staniland, P. Bucktrout, and A. G. Wood. 2023a. “Antarctic Fur Seal Abundance at South Georgia [Dataset].” https://doi.org/https://doi.org/10.5061/dryad.2rbnzs7tc.
———. 2023b. “Ninety Years of Change, from Commercial Extinction to Recovery, Range Expansion and Decline for Antarctic Fur Seals at South Georgia.” Global Change Biology 29 (24): 6867–87. https://doi.org/https://doi.org/10.1111/gcb.16947.
Forcada, J., and I. J. Staniland. 2018. “Antarctic Fur Seal: Arctocephalus Gazella.” In Encyclopedia of Marine Mammals (Third Edition), edited by B. Würsig, J. G. M. Thewissen, and K. M. Kovacs, Third Edition, 25–27. Academic Press. https://doi.org/https://doi.org/10.1016/B978-0-12-804327-1.00046-7.
Forcada, J., P. N. Trathan, and E. J. Murphy. 2008. “Life History Buffering in Antarctic Mammals and Birds Against Changing Patterns of Climate and Environmental Variation.” Global Change Biology 14 (11): 2473–88. https://doi.org/https://doi.org/10.1111/j.1365-2486.2008.01678.x.
Forcada, J., P. N. Trathan, K. Reid, and E. J. Murphy. 2005. “The Effects of Global Climate Variability in Pup Production of Antarctic Fur Seals.” Ecology 86 (9): 2408–17. https://doi.org/https://doi.org/10.1890/04-1153.
Garcia-Garin, O., I. García-Cuevas, M. Drago, D. Rita, M. Parga, M. Gazo, and L. Cardona. 2020. “No Evidence of Microplastics in Antarctic Fur Seal Scats from a Hotspot of Human Activity in Western Antarctica.” Science of The Total Environment 737: 140210. https://doi.org/https://doi.org/10.1016/j.scitotenv.2020.140210.
Gelman, A., Y. Goegebeur, F. Tuerlinckx, and I. van Van Mechelen. 2002. “Diagnostic Checks for Discrete Data Regression Models Using Posterior Predictive Simulations.” Journal of the Royal Statistical Society Series C: Applied Statistics 49 (2): 247–68. https://doi.org/10.1111/1467-9876.00190.
Gelman, A., and J. Hill. 2006. Data Analysis Using Regression and Multilevel/Hierarchical Models. Analytical Methods for Social Research. Cambridge University Press. https://doi.org/https://doi.org/10.1017/CBO9780511790942.
Gelman, A., J. Hwang, and A. Vehtari. 2014. “Understanding Predictive Information Criteria for Bayesian Models.” Statistics and Computing 24 (6): 997–1016. https://doi.org/10.1007/s11222-013-9416-2.
Gelman, A., and K. Shirley. 2011. “Inference from Simulations and Monitoring Convergence.” In Handbook of Markov Chain Monte Carlo, 162–74. https://doi.org/10.1201/b10905-7.
Gimenez, O., J. -D. Lebreton, R. Choquet, and R. Pradel. 2018. “R2ucare: An r Package to Perform Goodness-of-Fit Tests for Capture–Recapture Models.” Methods in Ecology and Evolution 9 (7): 1749–54. https://doi.org/https://doi.org/10.1111/2041-210X.13014.
Gimenez, O., V. Rossi, R. Choquet, C. Dehais, B. Doris, H. Varella, J. -P. Vila, and R. Pradel. 2007. “State-Space Modelling of Data on Marked Individuals.” Ecological Modelling 206 (3): 431–38. https://doi.org/https://doi.org/10.1016/j.ecolmodel.2007.03.040.
Goebel, M., D. Krause, S. Freeman, R. Burner, C. Bonin, R. Mercado, A. Van Cise, and J. Gafney. 2009. “Pinniped Research at Cape Shirreff , Livingston Island, Antarctica, 2008/09. In: AMLR 2008/2009 Field Season Report: Objectives, Accomplishments and Tentative Conclusions. NOAA-TM-NMFS-SWFSC-445.” https://doi.org/10.13140/RG.2.2.19451.64808.
Goebel, M., V. I. Vallejos, W. Z. Trivelpiece, R. S. Holt, and J. Acevedo. 2003. “Antarctic Fur Seal Pup Production in the South Shetland Islands.” Book Section. In AMLR 2001/2002 Field Season Report, NOAA-Technical Memorandum-NMFS-SWFSC-350., edited by J. Lipsky, 139–52. La Jolla, CA.
Goldsworthy, S. D., J. McKenzie, B. Page, M. L. Lancaster, P. D. Shaughnessy, L. P. Wynen, S. A. Robinson, K. J. Peters, A. M. M. Baylis, and R. R. McIntosh. 2009. “Fur Seals at Macquarie Island: Post-Sealing Colonisation, Trends in Abundance and Hybridisation of Three Species.” Journal Article. Polar Biology 32 (10): 1473–86. https://doi.org/10.1007/s00300-009-0645-y.
Gouhier, T., A. Grinsted, and Simko V. 2024. R Package biwavelet: Conduct Univariate and Bivariate Wavelet Analyses. https://github.com/tgouhier/biwavelet.
Grosbois, V., M. P. Harris, T. Anker-Nilssen, R. H. McCleery, D. N. Shaw, B. J. T. Morgan, and O. Gimenez. 2009. “Modeling Survival at Multi-Population Scales Using Mark–Recapture Data.” Journal Article. Ecology 90 (10): 2922–32. https://doi.org/https://doi.org/10.1890/08-1657.1.
Guinet, C., C. Barbraud, H. Weimerskirch, and K. Delord. 2025. “Monitoring Program Supported Financially and Logistically by the French Polar Institute IPEV (Projectgram 109 Ornithoeco, PI c. Barbraud/h. Weimerskirch), the Zone Atelier Antarctique (ZATA, CNRS-INEE), Terres Australes Et Antarctiques Françaises. This Study Is Part of the Long-Term Studies in Ecology and Evolution (SEE-Life) Program of the CNRS, c. Barbraud. All Work Was Carried Out in Accordance with the IPEV Ethics Committee.” Dataset.
Hill, Simeon L., Tony Phillips, and Angus Atkinson. 2013. “Potential Climate Change Effects on the Habitat of Antarctic Krill in the Weddell Quadrant of the Southern Ocean.” Journal Article. PLOS ONE 8 (8): e72246. https://doi.org/10.1371/journal.pone.0072246.
Hoffman, J. I., and J. Forcada. 2012. “Extreme Natal Philopatry in Female Antarctic Fur Seals (Arctocephalus Gazella).” Mammalian Biology 77 (1): 71–73. https://doi.org/https://doi.org/10.1016/j.mambio.2011.09.002.
Hoffman, J. I., P. N. Trathan, and W. Amos. 2006. “Genetic Tagging Reveals Extreme Site Fidelity in Territorial Male Antarctic Fur Seals Arctocephalus Gazella.” Molecular Ecology 15 (12): 3841–47. https://doi.org/https://doi.org/10.1111/j.1365-294X.2006.03053.x.
Hofmeyr, G. J. G. 2016. “Arctocephalus Gazella. The IUCN Red List of Threatened Species 2016: E.T2058A66993062. Http://Dx.doi.org/10.2305/IUCN.UK.2016-1.RLTS.T2058A66993062.en.” Electronic Article. http://dx.doi.org/10.2305/IUCN.UK.2016-1.RLTS.T2058A66993062.en.
Hofmeyr, G. J. G., M. N. Bester, A. B. Makhado, and P. A. Pistorius. 2006. “Population Changes in Subantarctic and Antarctic Fur Seals at Marion Island : Research Article.” South African Journal of Wildlife Research 36 (1): 55–68. https://doi.org/10.10520/EJC117228.
Hofmeyr, G. J. G., Marthán N. Bester, Pierre A. Pistorius, Tambu W. Mulaudzi, P. J. Nico de Bruyn, Justice A. Ramunasi, Hendrik N. Tshithabane, Trevor McIntyre, and Phathu M. Radzilani. 2007. “Median Pupping Date, Pup Mortality and Sex Ratio of Fur Seals at Marion Island.” South African Journal of Wildlife Research 37 (1): 1–8. https://doi.org/10.3957/0379-4369-37.1.1.
Hofmeyr, G. J. G., B. A. Krafft, S. P. Kirkman, M. N. Bester, C. Lydersen, and K. M. Kovacs. 2005. “Population Changes of Antarctic Fur Seals at Nyrøysa, Bouvetøya.” Polar Biology 28 (9): 725–31. https://doi.org/10.1007/s00300-005-0732-7.
Holdgate, M. W. 1963. “Fur Seals in the South Sandwhich Islands.” Polar Record 11 (73): 474–75. https://doi.org/https://doi.org/10.1017/S0032247400053742.
Hollyman, Trathan, P. 2021. “1Counts of Antarctic Fur Seals at Maiviken, South Georgia, from 2008 to 2020 (Version 1.0)[data Set].” NERC EDS UK Polar Data Centre. https://doi.org/https://doi.org/10.5285/41bfd98c-1a36-4689-a578-ad584426f9d2.
Humble, E., K. K. Dasmahapatra, A. Martinez-Barrio, I. Gregório, J. Forcada, A. -C. Polikeit, S. D. Goldsworthy, et al. 2018. “RAD Sequencing and a Hybrid Antarctic Fur Seal Genome Assembly Reveal Rapidly Decaying Linkage Disequilibrium, Global Population Structure and Evidence for Inbreeding.” G3 Genes|Genomes|Genetics 8 (8): 2709–22. https://doi.org/10.1534/g3.118.200171.
IUCN. 2024. “Standards and Petitions Committee. Guidelines for Using the IUCN Red List Categories and Criteria. Version 16. Prepared by the Standards and Petitions Committee. Downloadable from Https://Www.iucnredlist.org/Documents/RedListGuidelines.pdf.” https://www.iucnredlist.org/documents/RedListGuidelines.pdf.
Johnstone, G. W. 1982. “Part Three - Scientific Reports. Zoology. Expedition to the Australian Territory of Heard Island and the MacDonald Islands, 1980. Division of National Mapping, Technical Report 31.” Report.
Jones, K. A., N. Ratcliffe, S. C. Votier, S. Lisovski, A. -S. Bonnet-Lebrun, and I. J. Staniland. 2021. “Sexual Segregation in Juvenile Antarctic Fur Seals.” Journal Article. Oecologia 197 (2): 339–52. https://doi.org/10.1007/s00442-021-04983-y.
Jones, O. R., P. Barks, I. M. Stott, T. D. James, S. C. Levin, W. K. Petry, P. Capdevila, et al. 2022. “Rcompadre and Rage - Two R Packages to Facilitate the Use of the COMPADRE and COMADRE Databases and Calculation of Life History Traits from Matrix Population Models.” Methods in Ecology and Evolution 13: 770–81. https://doi.org/10.1111/2041-210X.13792.
Kawaguchi, So, Angus Atkinson, Dominik Bahlburg, Kim S. Bernard, Emma L. Cavan, Martin J. Cox, Simeon L. Hill, Bettina Meyer, and Devi Veytia. 2024. “Climate Change Impacts on Antarctic Krill Behaviour and Population Dynamics.” Journal Article. Nature Reviews Earth & Environment 5 (1): 43–58. https://doi.org/10.1038/s43017-023-00504-y.
Kellner, K. 2024. jagsUI: A Wrapper Around ’Rjags’ to Streamline ’JAGS’ Analyses. https://CRAN.R-project.org/package=jagsUI.
Kéry, M., and M. Schaub. 2012a. “Chapter 10 - Estimation of Survival, Recruitment, and Population Size from Capture–Recapture Data Using the Jolly–Seber Model.” In Bayesian Population Analysis Using WinBUGS, edited by M. Kéry and M. Schaub, 315–46. Boston: Academic Press. https://doi.org/https://doi.org/10.1016/B978-0-12-387020-9.00010-9.
———. 2012b. “Chapter 9 - Estimation of Survival and Movement from Capture–Recapture Data Using Multistate Models.” In Bayesian Population Analysis Using WinBUGS, edited by M. Kéry and M. Schaub, 263–313. Boston: Academic Press. https://doi.org/https://doi.org/10.1016/B978-0-12-387020-9.00009-2.
Kightley, S. P. J., and J. R. Caldwell. 1982. “The First Record of a Fur Seal Birth on Signy Island, South Orkney Islands.” British Antarctic Survey Bulletin 51: 287–89. https://nora.nerc.ac.uk/id/eprint/524417.
Krause, D. J., C. A. Bonin, M. E. Goebel, Christian S. Reiss, and George M. Watters. 2022. “The Rapid Population Collapse of a Key Marine Predator in the Northern Antarctic Peninsula Endangers Genetic Diversity and Resilience to Climate Change.” Frontiers in Marine Science Volume 8 - 2021. https://doi.org/10.3389/fmars.2021.796488.
Krause, D. J., R. L. Brownell Jr, C. A. Bonin, S. M. Woodman, D. Shaftel, and G. M. Watters. 2024. “Evaluating Threats to South Shetland Antarctic Fur Seals Amidst Population Collapse.” Mammal Review 54 (1): 30–46. https://doi.org/https://doi.org/10.1111/mam.12327.
Krause, D. J., S. M Woodman, and M. E. Goebel. 2025. “South Shetland Antarctic Fur Seal Pup Census. Version 1.1. SCAR - AntOBIS. Samplingevent Dataset. Https://Ipt-Obis.gbif.us/Resource?r=usamlr-Ssafs-Pup-Census&v=1.1.” Dataset. https://ipt-obis.gbif.us/resource?r=usamlr-ssafs-pup-census&v=1.1.
Lahoz-Monfort, J. J., B. J. T. Morgan, M. P. Harris, S. Wanless, and S. N. Freeman. 2011. “A Capture–Recapture Model for Exploring Multi-Species Synchrony in Survival.” Journal Article. Methods in Ecology and Evolution 2 (1): 116–24. https://doi.org/https://doi.org/10.1111/j.2041-210X.2010.00050.x.
Larsen, G. D., M. A. Cimino, J. Dale, A. S. Friedlaender, M. A. Goerke, and D. W. Johnston. 2025. “Terrestrial Spatial Distribution and Summer Abundance of Antarctic Fur Seals () Near Palmer Station, Antarctica, from Drone Surveys.” Ecology and Evolution 15 (4): e70833. https://doi.org/https://doi.org/10.1002/ece3.70833.
Laws, R. M. 1973. “Population Increase of Fur Seals at South Georgia.” Polar Record 16 (105): 856–58. https://doi.org/10.1017/S003224740006397X.
———. 1981. “Seal Surveys, South Orkney Islands, 1971 and 1974.” Journal Article. British Antarctic Survey Bulletin 54: 136–39. https://nora.nerc.ac.uk/id/eprint/524717.
Lea, M., C. Guinet, Y. Cherel, M. D. Hindell, L. Dubroca, and S. Thalmann. 2008. “Colony-Based Foraging Segregation by Antarctic Fur Seals at the Kerguelen Archipelago.” Journal Article. Marine Ecology Progress Series 358: 273–87. https://www.int-res.com/abstracts/meps/v358/p273-287/.
Lebreton, J. -D., and O. Gimenez. 2013. “Detecting and Estimating Density Dependence in Wildlife Populations.” The Journal of Wildlife Management 77 (1): 12–23. https://doi.org/https://doi.org/10.1002/jwmg.425.
Link, W. A., and R. J. Barker. 2005. “Modeling Association Among Demographic Parameters in Analysis of Open Population Capture–Recapture Data.” Biometrics 61 (1): 46–54. https://doi.org/10.1111/j.0006-341X.2005.030906.x.
Lowther, A. D., I. Staniland, C. Lydersen, and K. M. Kovacs. 2020. “Male Antarctic Fur Seals: Neglected Food Competitors of Bioindicator Species in the Context of an Increasing Antarctic Krill Fishery.” Scientific Reports 10 (1): 18436. https://doi.org/10.1038/s41598-020-75148-9.
Lunn, N. J., and I. L. Boyd. 1991. “Pupping-Site Fidelity of Antarctic Fur Seals at Bird Island, South Georgia.” Journal of Mammalogy 72 (1): 202–6. https://doi.org/ https://doi.org/10.2307/1381999.
March, D., M. Drago, M. Gazo, M. Parga, D. Rita, and L. Cardona. 2021. “Winter Distribution of Juvenile and Sub-Adult Male Antarctic Fur Seals (Arctocephalus Gazella) Along the Western Antarctic Peninsula.” Journal Article. Scientific Reports 11 (1): 22234. https://doi.org/10.1038/s41598-021-01700-w.
McBride, M. M., O. Schram Stokke, A. H. H. Renner, B. A. Krafft, O. A. Bergstad, M. Biuw, A. D. Lowther, and J. E. Stiansen. 2021. “Antarctic Krill Euphausia Superba: Spatial Distribution, Abundance, and Management of Fisheries in a Changing Climate.” Journal Article. Marine Ecology Progress Series 668: 185–214. https://www.int-res.com/abstracts/meps/v668/p185-214/.
McDonald, J. L., T. Bailey, R. J. Delahay, R. A. McDonald, G. C. Smith, and D. J. Hodgson. 2016. “Demographic Buffering and Compensatory Recruitment Promotes the Persistence of Disease in a Wildlife Population.” Ecology Letters 19 (4): 443–49. https://doi.org/https://doi.org/10.1111/ele.12578.
Meyer, W., B. G. Walker, and R. Holt. 1996. “Antarctic Fur Seal Abundance and Distribution in the South Shetland Islands, 1996.” Book Section. In AMLR 1995/96 Field Season Report (Vol. NOAA-AR-LJ 96-15). Antarctic Ecosystem Research Group. SWFSC/NMFS/NOAA., edited by J. Martin, NOAA-AR-LJ 96-15:97–99. SWFSC/NMFS/NOAA, La Jolla, California. https://library.oarcloud.noaa.gov/noaa_documents.lib/NMFS/SWFSC/AMLR_field_season_report/1995-1996.pdf.
Murphy, E. J., P. N. Trathan, J. L. Watkins, K. Reid, M. P. Meredith, J. Forcada, S. E. Thorpe, N. M. Johnston, and P. Rothery. 2007. “Climatically Driven Fluctuations in Southern Ocean Ecosystems.” Journal Article. Proceedings of the Royal Society B-Biological Sciences 274 (1629): 3057–67. https://doi.org/10.1098/rspb.2007.1180.
Ouled-Cheikh, J., D. March, R. Borras-Chavez, M. Drago, M. E. Goebel, J. M. Fariña, M. Gazo, M. Coll, and L. Cardona. 2024. “Future Climate-Induced Distribution Shifts in a Sexually Dimorphic Key Predator of the Southern Ocean.” Global Change Biology 30 (3): e17191. https://doi.org/https://doi.org/10.1111/gcb.17191.
Page, B., A. Welling, M. Chambellant, S. D. Goldsworthy, T. Dorr, and R. van Veen. 2003. “Population Status and Breeding Season Chronology of Heard Island Fur Seals.” Journal Article. Polar Biology 26 (4): 219–24. https://doi.org/10.1007/s00300-003-0478-z.
Paijmans, A. J., M. A. Stoffel, M. N. Bester, A. C. Cleary, P. J. N. De Bruyn, J. Forcada, M. E. Goebel, et al. 2020. “The Genetic Legacy of Extreme Exploitation in a Polar Vertebrate.” Scientific Reports 10 (1): 5089. https://doi.org/10.1038/s41598-020-61560-8.
Parmelee, D, W Fraser, BND Glass, and DR Neilson. 1977. “Ecological and Behavioral Adaptations to Antarctic Environments.” Journal Article. Antarct. JUS 12: 17.
Payne, M. R. 1977. “Growth of a Fur Seal Population.” Journal Article. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences 279 (963): 67–79. https://doi.org/DOI 10.1098/rstb.1977.0072.
———. 1979. “Growth in the Antarctic Fur Seal Arctocephalus-Gazella.” Journal Article. Journal of Zoology 187 (Jan): 1–20. <Go to ISI>://WOS:A1979GF56900001.
Plummer, M. 2003. “JAGS: A Program for Analysis of Bayesian Graphical Models Using Gibbs Sampling.” 3rd International Workshop on Distributed Statistical Computing (DSC 2003); Vienna, Austria 124 (April): 1–10.
———. 2024. Rjags: Bayesian Graphical Models Using MCMC. https://CRAN.R-project.org/package=rjags.
Pradel, R., and J. -D. Lebreton. 1999. “Comparison of Different Approaches to the Study of Local Recruitment of Breeders.” Bird Study 46 (sup1): S74–81. https://doi.org/10.1080/00063659909477234.
Pradel, R., C. M. A. Wintrebert, and O. Gimenez. 2003. “A Proposal for a Goodness-of-Fit Test to the Arnason-Schwarz Multisite Capture-Recapture Model.” Biometrics 59 (1): 43–53. https://doi.org/https://doi.org/10.1111/1541-0420.00006.
R_Core_Team. 2024. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.
Riecke, T. V., P. J. Williams, T. L. Behnke, D. Gibson, A. G. Leach, B. S. Sedinger, P. A. Street, and J. S. Sedinger. 2019. “Integrated Population Models: Model Assumptions and Inference.” Methods in Ecology and Evolution 10 (7): 1072–82. https://doi.org/https://doi.org/10.1111/2041-210X.13195.
Royle, J. A., and R. M. Dorazio. 2009. “10 - Modeling Population Dynamics.” In Hierarchical Modeling and Inference in Ecology, edited by J. A. Royle and R. M. Dorazio, 325–45. San Diego: Academic Press. https://doi.org/https://doi.org/10.1016/B978-0-12-374097-7.00012-0.
Royle, J. A., R. M. Dorazio, and W. A. Link. 2007. “Analysis of Multinomial Models with Unknown Index Using Data Augmentation.” Journal of Computational and Graphical Statistics 16 (1): 67–85. https://doi.org/10.1198/106186007X181425.
Santin-Janin, H., B. Hugueny, P. Aubry, D. Fouchet, O. Gimenez, and D. Pontier. 2014. “Accounting for Sampling Error When Inferring Population Synchrony from Time-Series Data: A Bayesian State-Space Modelling Approach with Applications.” Electronic Book Section. Public Library of Science. https://doi.org/10.1371/journal.pone.0087084.
Savoca, M. S., M. Kumar, Z. Sylvester, M. F. Czapanskiy, B. Meyer, J. A. Goldbogen, and C. M. Brooks. 2024. “Whale Recovery and the Emerging Human-Wildlife Conflict over Antarctic Krill.” Journal Article. Nat Commun 15 (1): 7708. https://doi.org/10.1038/s41467-024-51954-x.
SCAR. 2006. “Proposal to de-List Antarctic Fur Seals as Specially Protected Species. Working Paper 39. Antarctic Treaty Consultative Meeting XXIX, Edinburgh, UK, 12–23 June, 2006.” Report. https://scar.org/~documents/policy/antarctic-treaty/atcm-xxix-and-cep-ix-2006/atcm29-wp039?layout=default.
Schaub, M., O. Gimenez, B. R. Schmidt, and R. Pradel. 2004. “Estimating Survival and Temporary Emigration in the Multistate Capture-Recapture Framework.” Ecology 85 (8): 2107–13. https://doi.org/https://doi.org/10.1890/03-3110.
Schwarz, L. K., M. E. Goebel, D. P. Costa, and A. M. Kilpatrick. 2013. “Top-down and Bottom-up Influences on Demographic Rates of Antarctic Fur Seals Arctocephalus Gazella.” Journal Article. Journal of Animal Ecology 82 (4): 903–11. https://doi.org/10.1111/1365-2656.12059.
Siniff, Donald B., Robert A. Garrott, Jay J. Rotella, William R. Fraser, and David G. Ainley. 2008. “Opinion: Projecting the Effects of Environmental Change on Antarctic Seals.” Journal Article. Antarctic Science 20 (05): 425–35. https://doi.org/doi:10.1017/S0954102008001351.
Vehtari, A., A. Gelman, D. Simpson, B. Carpenter, and P. -C. Bürkner. 2021. Rank-Normalization, Folding, and Localization: An Improved \(\widehat{R}\) for Assessing Convergence of MCMC (with Discussion).” Bayesian Analysis 16 (2): 667–718. https://doi.org/10.1214/20-BA1221.
Waluda, C. M., S. Gregory, and M. J Dunn. 2010. “Long-Term Variability in the Abundance of Antarctic Fur Seals Arctocephalus Gazella at Signy Island, South Orkneys.” Polar Biol. 33 (3): 305–12.
Watanabe, Sumio. 2010. “Asymptotic Equivalence of Bayes Cross Validation and Widely Applicable Information Criterion in Singular Learning Theory.” Journal of Machine Learning Research 11 (116): 3571–94. http://jmlr.org/papers/v11/watanabe10a.html.
Watters, George M., Jefferson T. Hinke, and Christian S. Reiss. 2020. “Long-Term Observations from Antarctica Demonstrate That Mismatched Scales of Fisheries Management and Predator-Prey Interaction Lead to Erroneous Conclusions about Precaution.” Journal Article. Scientific Reports 10 (1): 2314. https://doi.org/10.1038/s41598-020-59223-9.
Weegman, M. D., T. W. Arnold, R. G. Clark, and M. Schaub. 2021. “Partial and Complete Dependency Among Data Sets Has Minimal Consequence on Estimates from Integrated Population Models.” Ecological Applications 31 (3): e2258. https://doi.org/https://doi.org/10.1002/eap.2258.
Wege, M., M. -P. Etienne, C. W. Oosthuizen, R. R. Reisinger, M. N. Bester, and P. J. N. de Bruyn. 2016. “Trend Changes in Sympatric Subantarctic and Antarctic Fur Seal Pup Populations at Marion Island, Southern Ocean.” Marine Mammal Science 32 (3): 960–82. https://doi.org/https://doi.org/10.1111/mms.12306.
Wickens, P., and A. E. York. 1997. “Comparative Population Dynamics of Fur Seals.” Marine Mammal Science 13 (2): 241–92. https://doi.org/https://doi.org/10.1111/j.1748-7692.1997.tb00631.x.
Wood, S. N. 2016. “Just Another Gibbs Additive Modeler: Interfacing JAGS and Mgcv.” Journal of Statistical Software 75 (7): 1–15. https://doi.org/10.18637/jss.v075.i07.